Appendix 3. 



EXPERIMENTAL STUDffiS 

OF A SIDE IMPACT INFLATABLE HEAD PROTECTOR 

2. INTRODUCTION 

This paper describes experimental studies of a side inflatable head protector used in 
"Satum"cars. The main objective of the studies is to find out relationships between the pressure 
within the airbag, its elongation and the tension force in tethers attached to the car's body that 
extend the airbag longitudinally. 

3. OBJECT OF STUDIES 

The side protection airbag for "Saturn" (Fig. 73) is an elongated construction that consists of two 
flat sheets of a strong fabric impregnated with an elastic synthetic material and glued along their 
perimeters and in interior areas. The figure is a photogrammetry: sizes can be measured on it to 
scale. Four reference points on the figure are placed in the vertices of a rectangle 60 cm by 80 
cm. 




Fig. 73. A photogrammetry of the side head protector for "Saturn" 



The airbag is equipped with two fixing tethers. The airbag is placed along the side wall of the car 
body and is attached to the body along its upper edge above the car door openings; also, it is at- 
tached to the front and back parts of the body with the tethers. The tethers are needed to create a 
tension in the airbag and thus reduce its deformability. 

4. GOAL OF STUDIES 

As the airbag is inflated, cylindrical convexities tilted to its longitudinal axis appear in its 
working area (Fig. 74). The length of the bag is thus reduced. As the tethers resist to it, a tension 
force appears in the airbag. 

The tension force depends on the pressure within the airbag and on the distance between points 
where the tethers are attached. The goal of this study is to find out these dependences experimen- 
tally. 
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Fig. 74. The appearance of the airbag when inflated. 

5. 3KCIIEPHMEHTAJIbHAflyCTAHOBKA ® 

The experimental bench (Fig. 75) includes a Base, a Tensioning device, an Inflating device, a 
Dy-namometer and a Pressure gauge. 




'ig. 75. The 



general view of the experimental bench. (1) the base, (2) the tensioning device, 
(3) the inflating device, (4) the dynamometer, (5) the pressure gauge 

The base (Fig. 76) is made of a wooden plank 50 mm thick. Blocks of wood are attached to the 
plank using screws and glue to make a support for the tensioning device and the dynamometer. 
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Fig. 76. The base of the experimental bench 
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The tensioning device (Fig. 77) is a bar with a rack along which moves a carriage with a pinion. 

The pinion is toothed with the rack and has a handle. The position of the carriage on the rack can be 
measured using a centimeter scale on the bar and fixed using a lock screw. The tensioning device is 
actually a post from a photographic enlarger "Belarus". 




Fig. 77. The tensioning device 

The inflating device (Fig. 78) is a pump with a valve and a pressure gauge to control the inflation. 
The inflating device is actually a standard foot-operated car pump. 




Fig. 78. The inflating device 



The dynamometer (Fig. 79) is a pointer instrument for force measurements graduated for kilograms 
of force with the measurement Umits 0 to 20 kgf. Its face has the diameter of about 20 cm. The 
dynamometer has an ear to hang a weight or to apply a force. 
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Fig. 79. The dynamometer 

The pressure gauge (Fig. 80) is a pointer device for pressure measurements graduated for kgf^cm^ 
with the measurement Umits 0 to 1.6 kgfi^cm^ . The diameter of its face is about 16 cm. 



One tether of the airbag is attached to the dynamometer's hook and the other one to the carriage of 
the tensioning device using a steel stripe with ears. 

A bronze cork with a socket fabricated specially for this purpose is inserted into the air hole of the 
airbag. A rubber hose is connected to the socket. The other end of the hose is connected to a tee 
joint. The pressure gauge and the inflating device are connected to the same T-joint. 

6. DETERMINING THE ELASTIC COMPLIANCE OF THE DYNAMOMETER 

Components of the experimental bench are stiff, except for the dynamometer. Their compliance, 
comparing to that of the airbag, can be neglected a priori. The only doubt was about the 
dynamometer that is noticeably compliant. To measure this comphance quantitatively, two series of 
measurements were performed. The traction of the dynamometer and die displacement of its ear 
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were measured. The displacements were measured by a clock-type indicator (Fig. 81) with the 
measurement limits 0 to 1 mm and the scale factor 0,001 mm. 




Fig. 81. A clock-type indicator 

The table of the measurement results is given below. 

Tab. 1. Numerical results of the dynamometer compliance investigation 



First series 


Second series 


Fkgf 


S mm 


Fkgf 


S mm 


Fkgf 


S mm 


6.0 


0.058 


2.0 


0.049 


5.0 


0.354 


6.2 


0.080 


2.2 


0.068 


6.0 


0.450 


6.4 


0.100 


2.4 


0.090 


7.0 


0.550 


6.6 


0.122 


2.6 


0.111 


8.0 


0.660 


6.8 


0.141 


2.8 


0.131 






7.0 


0.160 


3.0 


0.150 






7.2 


0.184 


3.2 


0.172 






7.4 


0.202 


3.4 


0.192 






7.6 


0.224 


3.6 


0.212 






8.0 


0.269 


3.8 


0.230 






8.2 


0.287 


4.0 


0.250 






8.4 


0.305 


4.2 


0.272 






8.6 


0.329 


4.4 


0.290 






8.8 


0.349 


4.6 


0.310 






9.0 


0.372 


4.8 


0.330 
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On the basis of these numerical data, hnear dependencies between the dynamometer traction and the 
position of its moving ear were built by the method of least squares. Those are shown as solid lines 
on Fig. 82. 
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Fig. 82. Results of the dynamometer investigation 

The compliance factor for the first series of measurements was 0.104 mm/kgf, and 0.101 mm/kgf for 
the second series. The additional displacement caused by the dynamometer compliance during 
measurements of the airbag with the traction of up to 19 kgf does not exceed 2 nmi that is about 1 to 
2% of the airbag extension. Therefore we can neglect the influence of the dynamometer compliance 
on the measurement results. 

7. AIRBAG TESTING TECHNIQUE 

The technique of testing was as follows. The airbag was fixed on the experimental bench so that the 
distance between the tether fixation points was a predefined small number. Then the pressure in the 
airbag was maximized. The tension force in the airbag was measured. Then the pressure in the 
airbag was reduced to the next fixed value. Again, the tension force was measured. 

Notice that the airbag is not hermetical, and the pressure within falls naturally. The venting places 
are located in concavities of the large hollows (the A area on Fig. 73). According to rough estimates, 
the permeability is such that the airbag permits about 0.1 liter per second to escape at the pressure of 
0.6kg^cm^ 

Making the pressure in the airbag lower step by step and measuring the airbag tension at each step, 
we were obtaining a curve of dependence between the tension force and the internal pres-sure for a 
fixed distance between the tether fixation points. 

After that, the distance between the tether fixation points was increased by 1 to 2 cm, and an-other 
series of measurements was done. Such series followed one another until the distance between the 
fixation points became maximal. 

Full sets of measurements were repeated many times in different days. 

Before each series the length of the empty airbag lo was measured. The length of the empty air-bag 
means the distance between the centers of fixation holes of the aiibag tethers measured at the tension 
force of 0.5 kgf. The empty airbag length increased noticeably after the first series of tests. After that 
it changed little. 




Page 97 



8, RESULTS OF AIRBAG TESTS 

Below, tables of measurement results are given. Fields of the tables contain the registration of the 
dynamometer that measured the tension force /in kgf. At the head of each column the regis-tration 
of the pressure gauge stands: pressure above ambient p in the airbag in kgfi^cm^ . The left cell of each 
row contains the registration of the tensioning device scale x in cm. 



Tab. 2. Tension/in kgf. First series of measurements 1 1/1 1/01. 
Empty bag length before measurement /o= 183.9 cm 





p, kgfi'cm^ 


xcm 


0.02 


0.04 


0.06 


0.08 


0.10 


0.12 


0.14 


0.16 


0.18 


0.20 


2 


0.00 


0.10 


0.20 


0.25 


0.30 


0.38 


0.42 


0.50 


0.58 


0.60 


4 


0.05 


0.17 


0.23 


0.33 


0.42 


0.55 


0.65 


0.72 


0.90 


0.95 


6 


0.12 


0.22 


0.33 


0.43 


0.60 


0.70 


0.81 


0.95 


1.03 


1.18 


8 


0.15 


0.30 


0.45 


0.62 


0.93 


1.00 


1.15 


1.30 


1.42 


1.58 


10 


0.20 


0.35 


0.43 


0.65 


0.85 


1.05 


1.21 


1.39 


1.58 


1.78 


12 


0.20 


0.40 


0.88 


1.10 


1.40 


1.70 


2.00 


2.23 


2.50 


2.70 


14 


0.25 


0.50 


1.00 


1.45 


1.95 


2.33 


2.65 


3.00 


3.40 


3.70 


15 


0.30 


0.65 


1.30 


1.90 


2.45 


2.90 


3.30 


3.80 


4.20 


4.50 


16 


0.35 


0.90 


1.65 


2.35 


2.93 


3.50 


4.10 


4.60 


5.10 


5.50 


17 


0.40 


1.15 


2.05 


2.90 


3.65 


4.26 


4.90 


5.55 


6.10 


6.50 


18 


0.50 


1.45 


2.57 


3.50 


4.35 


5.00 


5.75 


6.40 


6.95 


7.60 


19 


0.62 


1.80 


3.00 


4.10 


5.10 


5.87 


6.70 


7.60 


8.20 


8.80 


20 


0.70 


2.10 


3.50 


4.80 


5.80 


6.85 


7.80 


8.60 


9.40 


10.20 


21 


0.90 


2.58 


4.20 


5.50 


6.70 


8.00 


8.80 


9.70 


10.80 


11.60 


22 


1.00 


2.90 


4.70 


6.30 


7.40 


8.80 


10.00 


11.10 


12.10 


13.10 


23 


1.35 


3.40 


5.40 


7.10 


8.60 


9.90 


11.10 


12.50 


13.70 


14.70 


24 


1.70 


3.95 


6.20 


8.00 


9.60 


11.20 


12.60 


14.00 


15.40 


16.60 



Tab. 3. Tension/in kgf. Second series of measurements 11/13/01 
Empty bag length before measurement /o=186.8 cm 





p, kgfcm^ 


X cm 


0.02 


0.04 


0.06 


0.08 


0.10 


0.12 


0.14 


0.16 


0.18 


0.20 


10 


0.20 


0.40 


0.60 


0.80 


1.10 


1.25 


1.40 


1.60 


1.80 


2.00 


12 


0.30 


0.60 


1.00 


1.35 


1.65 


1.95 


2.20 


2.50 


2.72 


3.00 


14 


0.40 


0.80 


1.30 


1.90 


2.30 


2.70 


3.00 


3.40 


3.75 


4.10 


16 


0.50 


1.20 


1.90 


2.60 


3.20 


3.70 


4.30 


5.00 


5.50 


6.00 


18 


0.85 


1.60 


2.86 


3.70 


4.60 


5.40 


6.20 


6.70 


7.50 


8.00 
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20 


1.22 


2.70 


4.00 


5.20 


6.20 


7.40 


8.30 


9.30 


10.10 


11.00 


22 


1.80 


3.60 


5.30 


6.80 


8.30 


9.50 


11.00 


12.00 


13.30 


14.30 


24 


2.60 


5.00 


7.00 


9.00 


11.00 


12.40 


14.20 


15.70 


17.20 


18.40 


Tab. 4. Tension/in kgf. Third series of measuiements 11/15/01 
Empty bag length before measurement iO= 186.8 cm 






gf^cm^ 


X cm 


0.02 


0.04 


0.06 


0.08 


0.10 


0.12 


0.14 


0.16 


0.18 


0.20 


10 


0.10 


0.25 


0.50 


0.70 


0.90 


1.05 


1.25 


1.40 


1.62 


1.80 


12 


0.20 


0.41 


0.75 


1.05 


1.40 


1.65 


2.00 


2.25 


2.50 


2.70 


14 


0.28 


0.65 


1.15 


1.60 


2.15 


2.55 


2.90 


3.30 


3.65 


4.00 


16 


0.40 


1.10 


1.85 


2.50 


3.20 


3.65 


4.22 


4.80 


5.30 


5.70 


18 


0.70 


1.70 


2.70 


3.60 


4.30 


5.20 


6.00 


6.70 


7.30 


7.90 


20 


1.00 


2.40 


3.70 


5.00 


6.00 


7.10 


8.10 


9.00 


9.90 


10.80 


22 


1.50 


3.30 


5.10 


6.70 


8.00 


9.40 


10.60 


12.00 


13.00 


14.00 


24 


2.10 


4.50 


6.80 


8.70 


10.40 


12.20 


13.90 


15.40 


16.70 


18.40 



Tab. 5. Tension/in kgf. Fourth series of measurements 11/20/01 
Empty bag length before measurement iO= 186.8 cm 





p, kgfi'cm^ 


X cm 


0.02 


0.04 


0.06 


0.08 


0.10 


0.12 


0.14 


0.16 


0.18 


0.20 


10 


0.19 


0.43 


0.70 


0.95 


1.20 


1.50 


1.60 


1.90 


2.00 


2.1 


12 


0.30 


0.70 


1.00 


1.30 


1.65 


2.00 


2.20 


2.50 


2.70 


3.10 


14 


0.40 


0.90 


1.50 


2.00 


2.36 


2.80 


3.20 


3.60 


4.00 


4.30 


16 


0.50 


1.20 


1.90 


2.50 


3.10 


3.70 


4.40 


5.00 


5.40 


5.80 


18 


0.70 


1.70 


2.70 


3.60 


4.50 


5.30 


6.00 


6.80 


7.50 


8.00 


20 


1.00 


2.40 


3.90 


5.10 


6.20 


7.10 


8.20 


9.10 


10.00 


10.90 


22 


1.50 


3.40 


5.30 


6.80 


8.20 


9.50 


10.90 


11.90 


13.10 


14.20 


24 


2.00 


4.60 


7.00 


8.90 


10.70 


12.40 


13.90 


15.50 


17.00 


18.30 



Tab. 6. Tension/in kgf. Fifth series of measurements 11/20/01 
Empty bag length before measurement iO=186.8 cm 







gf^cm^ 


X cm 


0.02 


0.04 


0.10 


0.12 


0.16 


0.20 


0.30 


0.40 


0.50 


0.60 


12 


0.20 


0.50 


1.40 






2.70 


3.6 


4.30 




5.90 


14 


0.40 


0.80 


2.10 






3.80 


5.30 


6.80 


7.70 


8.50 
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16 


0.60 


1.20 




3.60 


0.45 


5.30 


7.50 


9.70 


11.00 


12.00 


18 


0.80 


1.70 


4.40 






9.90 


10.70 


13.00 


14.60 




Tab. 7. Tension/in kgf. Sixth series of measurements 1 1/21/01 
Empty bag length before measurement /o=186.8 cm 




p, kgf^cm^ 


jccm 


0.02 


0.04 


0.06 


0.08 


0.10 


0.12 


0.14 


0.16 


0.18 


0.20 


10 


0.15 


0.35 


0.60 


0.75 


1.00 


1.20 


1.30 


1.50 


1.70 


1.80 


12 


0.20 


0.60 


0.95 


1.30 


1.50 


1.90 


2.20 


2.40 


2.60 


2.80 


14 


0.35 


0.80 


1.40 


1.80 


2.10 


2.60 


3.10 


3.50 


3.80 


4.10 


16 


0.40 


1.20 


1.80 


2.50 


3.10 


3.70 


4.35 


4.90 


5.30 


5.80 


18 


0.60 


1.70 


2.50 


3.60 


4.30 


5.20 


6.20 


6.90 


7.60 


8.00 


20 


0.90 


2.40 


3.80 


4.90 


6.10 


7.20 


8.10 


9.10 


10.00 


10.80 


22 


1.40 


3.30 


5.30 


6.80 


8.20 


9.40 


10.90 


12.10 


13.10 


14.00 


24 


2.00 


4.40 


6.9 


8.90 


10.60 


12.30 


13.90 


15.50 


17.00 


18.30 



Tab. 8. Tension/in kgf. Seventh series of measurements at a higher pressure 11/22/01 

Empty bag length before measurement /o=186.8 cm 





p, kgfi'cm^ 


xcm 


0.1 


0.2 


0.3 


0.4 


0.5 


0.6 


10 


1.1 


2.0 


2.8 


3.5 


4.1 


4.6 


12 


1.4 


2.6 


3.6 


4.6 


5.4 


6.0 


14 


2.1 


3.9 


5.4 


6.7 


7.7 


8.6 


16 


3.1 


5.6 


7.5 


9.3 


10.7 


12.0 


18 


4.2 


7.6 


10.2 


12.5 


14.5 


16.2 


20 


5.6 


9.8 


13.2 


16.0 


18.7 


20.6 



Results of the first series are sUghtly different from those of subsequent series. This difference was 
apparently caused by non-elastic deformations of the airbag elements during the first series. Also, 
the testing technique was still under development in the course of the first series of measurements. It 
was noticed that the measurement of small forces can be substantially distorted due to the friction of 
the heavy cork and the airbag at the laboratory table. After that we decided to hang up the cork and 
stir the airbag from time to time to reduce the fiiction. The first series of measurements lasted for 
several hours, dierefore in this case the influence of creep of the tether and the airbag material was 
stronger. 

Nevertheless, actual readings of all measurement instruments are included in this report to illustrate 
the results and to facilitate the comparison. 
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9. MEASUREMENT RESULTS PROCESSED 

Two sets of curves have been built after results of the studies. The first set was built on the basis of 
measurements described in Tab. 3, Tab. 4, Tab. 5 and Tab. 7, and the second set on the basis of 
those described in Tab. 8, 

The sets of curves describe a functional dependence of the tension force in the airbag on the pressure 
above ambient within it and on the location of the tether attachment points. 

The following experimental readings of the measuring instruments were processed: 

/, - the tension of the airbag; 

p, the pressure above ambient widiin the airbag; 

X, the readings of the tensioning device scale that define the locations of the tether 
attachment points. 

The final diagrams use the argument d instead of the argument x: the difference of the distance 
between the points of the tether attachment to the experimental bench and the length of the empty 
airbag. This investigation refers to negative values of d (the distance between the tether attachment 
points is less than the empty airbag length). 

The distance between the tether attachment points in the experimental bench / and the tensioning 
device readings jc are related as follows: 

( 2 ) /= x+160.2 cm. 

Thus, 

( 3 ) d-x-k +160.2 cm; 

where lois the empty airbag lengdi before the test. 

Results of the first measurement series listed in Tab. 2 are shown on 

Fig. 83. 

Most representative data are found in Tab, 3, Tab. 4, Tab. 5 and Tab. 7. These tables contain data of 
measurements of the airbag tension in the range of the internal pressure above ambient p 0.02 to 0.2 
kgf^cm^ and with the x value 10 to 24 cm. In essence, those are results of the same measurement 
series performed at different times. The results correlate pretty well: it can be seen from diagrams of 
Fig. 84 that describe four series of experiments. 
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Fig. 83. Results of the first trial set of measurements 




Fig. 84. A comparison between the measurement results 
Grid colors: red for Tab. 3, green for Tab, 4, blue for Tab. 5, black for Tab. 7 7. 

An approximation of the dependence by a surface of the fourth order gives results shown in 
Tab, 9, 

Tab. 9. Averaged data firom Tab. 3, Tab. 4, Tab. 5, Tab. 7. 

/?, kgfi^cm^ 
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xcm 


0.02 


0.04 


0.06 


0.08 


0.10 


0.12 


0.14 


0.16 


0.18 


0.20 


10 


0.179 


0.383 


0.586 


0.795 


1.011 


1.231 


1.446 


1.640 


1.794 


1.882 


12 


0.247 


0.589 


0.910 


1.223 


1.532 


1.838 


2.137 


2.419 


2.667 


2.861 


14 


0.305 


0.837 


1.325 


1.784 


2.225 


2.653 


3.067 


3.461 


3.826 


4.143 


16 


0.424 


1.203 


1.910 


2.564 


3.180 


3.768 


4.332 


4.871 


5.377 


5.840 


18 


0.654 


1.741 


2.724 


3.625 


4.465 


5.258 


6.012 


6.729 


7.409 


8.043 


20 


1.026 


2.487 


3.807 


5.013 


6.129 


7.175 


8.162 


9.098 


9.986 


10.823 


22 


1.549 


3.455 


5.178 


6.750 


8.200 


9.551 


10.820 


12.019 


13.155 


14.230 


24 


2.215 


4.640 


6.835 


8.839 


10.685 


12.399 


14.003 


15.514 


16.943 


18.295 


y changing to variables p, <i,/we obtain Tab. 10. 

Tab. 10. Averaged data in variables p, d,f. 




p, kgf^cm^ 


d cm 


0.02 


0.04 


0.06 


0.08 


0.10 


0.12 


0.14 


0.16 


0.18 


0.20 


-14 


0.262 


0.655 


1.021 


1.373 


1.717 


2.056 


2.386 


2.699 


2.979 


3.208 


-12 


0.331 


0.932 


1.480 


1.992 


2.481 


2.952 


3.407 


3.842 


4.246 


4.604 


-10 


0.479 


1.344 


2.128 


2.850 


3.529 


4.173 


4.790 


5.378 


5.933 


6.444 


-8 


0.750 


1.942 


3.019 


4.006 


4.923 


5.786 


6.605 


7.384 


8.122 


8.813 


-6 


1.167 


2.754 


4.187 


5.497 


6.707 


7.838 


8.904 


9.915 


10.873 


11.777 


-4 


1.735 


3.789 


5.645 


7.340 


8.902 


10.355 


11.719 


13.007 


14.226 


15.380 


-2 


2.438 


5.034 


7.386 


9.533 


11.510 


13.344 


15.060 


16.675 


18.201 


19.645 


0 


3.239 


6.458 


9.381 


12.053 


14.512 


16.791 


18.918 


20.914 


22.797 


24.578 



The same information is shown visually as a 3D diagram on Fig. 85. 



Page 103 



Fig. 85. Averaged data in variables p, d^f, 
A set of curves shown on Fig. 86 is more convenient to use. 
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Fig. 86. A set of curves to determine the airbag tension force/ 

Results of the airbag test at an increased pressure are listed in Tab. 6 and Tab. 8. Systematic data are 
given in Tab. 8 that is used to obtain final results. Those are shown visually on Fig. 87. The same 
figure shows a surface that approximate the results. Tab. 1 1 lists ordinates of this surface. 
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Fig. 87. Results of measurements at the increased pressure and their approximation 

The surface that approximates the experimental data is shown in red. The difference is hardly no 
ticeable. Fig. 88 shows the same function as a set of curves. 



24 
22 
20 
18 
16 



0.1 



0.2 

d=-14 cm 
d=-12cm 
d=-10cm 
d= -8 cm 
d= -6 cm 



0.3 0.4 
pressure kgffsq.cm 



0.5 



0.6 



Fig. 88. A set of curves to determine the airbag tension force /at the increased pressure 
Tab. 11. Averaged data of measurement at the increased pressure 
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20 
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16.087 


18.646 


20.917 



It is noticeable that there is a substantial difference between the approximation functions on Fig. 86 
and Fig. 88. This difference was apparently caused by the airbag material creep that revealed itself 
when testing the airbag at different speeds. 

10. AIRBAG MATERIAL CREEP 



The following experiment was run to obtain tentative data about the material creep rate. The air-bag 
with no gas within was quickly stretched by 8 cm and maintained in this state until the stress relaxed 
(20 min). Then the airbag was brought back to its initial state. 

Tab. 12 that describes this experiment contains a row that lists reference time t in minutes from the 
beginning of the experiment, a row of displacements x in cm and a row of the tension force / 
measured in kilograms of force. 



Tab. 12. Results of the stress relaxation experiment 
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Fig. 89. A hysteresis loop observed in the airbag when loaded and unloaded 

The experiment shows that the tension in the airbag material relaxes by approximately the following 
value: 

( 4 ) (1 - 1 1.7/16.6) 100% = 30%. 
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The relaxation time is about 20 min. This explains, apparently, the discrepancy between results of 
experiments where the measurements were done at different speeds. 

11. CONCLUSIONS 

1. This report describes qualitative and quantitative results that characterize the dependence 
between the tension force, the internal pressure, and the distance between the tether fixation 
points (at relatively low pressures). 

2. It has been established that plastic deformations of the airbag material at the intemal pressure up 
to 0.6 kgf^cm2 and the tension up to 20 kgf have little effect on its elastic response. Though, 
deformations caused by creeping of the material affect the parameters of interest very 
essentially. 

3. The airbag is not hermetic and it permits the air to escape through seams where the flat sheets 
are connected. The venting places are located in concavities of the large hollows (the A area Fig. 
73). The permeability is 0.1 liter per second at the pressure of 0.6 kgfi^cm2. 

4. Designers of the airbag should try to avoid large hollows with concavities at the edge. 
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Side impact and roll-over airbags are similar to air mattresses attached to components of a car's 
body. The airbag consists essentially of two flat pieces of film glued between each other both along 
their contours and in intermediate areas. Areas where the sheets are connected have complicated 
shapes which can be seen, for instance, from the picture of a side impact airbag for the "Saturn" car 
shown on Fig. 11.1. 




Fig. 1 1 .1 . A side impact airbag for the "Saturn" car 



In convex sections of the glued areas, there is a concentration of stress that tears off the glued seam. 
A bigger radius of the rounded comers reduces this stress but at the same time decreases hollows 
filled with gas thus deteriorating the effectiveness of the airbag. So here we have a reason to 
formulate an optimization problem. The peeling strength of the adhesive seam may prove a 
determinant for the whole airbag design being an active limit for any optimizing solution. 

Before discussing an optimal design, one should learn how to determine a shape and stress of 
studied objects. 

The present work solves a number of simplified tasks where typical elements of sheets pasted 
together are contained. This problem is solved partly using analytical methods and partly 
numerically. The latter case is based on the finite elements method, and the computation was carried 
out using ANSYS5.5. 

12. AIRBAG WITH HOLLOWS PERPENDICULAR TO THE TENSION 

The peculiarity of the side impact airbag's design is the presence of a number of cylindrical hollows 
(pockets) filled with gas and arranged perpendicularly or at an angle to the longitudinal axis of the 
airbag. 

Seeing the complex configuration of the real side impact airbag, we will investigation simplified 
models that keep most relevant features of the real object of study. 

This part deals with a simplified model where the cylindrical pockets are arranged perpendicularly 
to the longitudinal direction of the airbag. The puipose of this study is to find relationships between 
the gas pressure, the airbag material tension and its elongation. 

The part contains differential equations derived, solutions thereof in the analytical form, and 
families of curves for a set of source data as close to the reality as possible. 
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12.1. Analytical model 

Let's suppose that the system contains an infinite sequence of cylindrical hollows with a constant 
spacing, and each hollow is infinitely long. The system is built of two infinite film pieces glued with 
stripes 2b wide. Between each couple of adjacent glued stripes there is a non-glued space 2a wide 
that forais a cylindrical hollow to be filled with gas under the pressure p. The thin film that the 
airbag is made of has the thickness h, elasticity modulus E and Poisson ratio v. The film is assumed 
to be absolutely flexible. 

Fig. 12.1 shows an excision from the infinite system that has the size one along the generatrix. This 
excision contains two cylindrical hollows and their adjacent glued stripes of the width b. Sizes a, b 
and size one refer to the non-deformed state of the film. When loaded, the film stretches and the 
sizes change. 



Fig. 12.1 An excision from a system with a sequence of infinite cylindrical hollows perpendicular to 

the longitudinal direction of the airbag 

As one can see from the above-stated, the shape of the bag can be drawn by a rectilinear generatrix 
moving along a curvilinear guide. 

The panel formed by the generatrices is in tension, and the linear distributed loading on it is 2q, 

The system is symmetric, so we can analyze only a repeatable elementary cell instead of the whole 
structure [ 19 ]. 



2( 
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A 




Fig. 12.2. An elementary cell for the system from Fig. 12.1 

Fig. 12.2 shows fragments of the film included in the elementary cell as bold lines. Note that only 
one of the glued films is included in the cell, and also only the quarter of the whole gas volume that 
fills one hollow. The elementary cell accepts only the half of the linear force that stretches the 
airbag. 

12.2. Mathematical formulation 

Let's place the elementary cell in a Cartesian coordinate system jci0;c2 as shown on Fig. 12.2. Then 
introduce the parameter t the physical meaning of which is the length of the guiding curve measured 
along the non-deformed surface of the film. The starting point of the measurement is the middle of 
the non-glued film area (Fig. 12.2). 

The variation limits of the parameter are: 

(12.1) 0<f<a+/?; 

where a is the half-width of the non-glued airbag area in the non-deformed state; 

h is the half-width of the glued airbag area in the non-deformed state. 

Value r= a corresponds to the boundary between the non-glued and glued areas. 

The airbag is a cylindrical body with a rectilinear generatrix, therefore its shape is completely 
defined by the firagment of the generatrix described above. 

Let this line be defined parametrically using coordinates jci(f), %%{{) that depend on the parameter f. 

The displacement method is used to solve the problem. New functions will be introduced during the 
derivation of the respective formulae. All of them will be functions of the parameter t. 

As all variables are ftmctions of one parameter r, it will be omitted below as an argument of 
functions. Differentiation with respect to this parameter will be denoted by a stroke. 

Obviously, the curve length element in the deformed state is: 
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(12.2) ds = dt4{x[f+{Af . 

The relative elongation of tiie film along tiie guide curve Ei can be described by the formula: 

(12.3) e^=^Z^ = /_i; 

where 

(12.4) s^4{x[f+{x[f . 

The relative elongation along tfie generatrix will be denoted as £2* Note that this is a constant tiiat 
does not depend on the parameter t. 

To determine the elastic potential of the system and find stresses in the film, we will need the 
physical relationships: 

E 

(12.5) =- ^-(ei+vej; 

1-v 

( 12.6) =]— ^(^^1+^2); 

where Gi and G2 are stresses along the generatrix and the guide, respectively. 
12.3. Full energy of the system 

We use the variation approach to solve the problem using the displacement method. This approach 
gives us a good formalization of this solution stage, helps choose a consistent system of variables, 
formulate natural boundary conditions, pose the dual problem if necessary. Also, there is no need to 
compose equilibrium equations looking at figures. 

Let's choose a system of functions and constants to completely define the deformed state of the 
system. These functions can be JCi, X2, and the constant £2. 

The full energy of the system W, that includes the energy of loading also, is now described by the 
following fimctional: 

Eh 



(12.7) W=\ 



0 L 



2(1 -v') 



[z[ + 2v£i£2 + £2 )- p%{a - t)xX (1 + £2 ) 



dt; 



where x(z) is the Heviside function equal to 1 at z > 0 and 0 at z < 0. 

The first term in the brackets is the elastic energy of the airbag material, and the second term is the 
energy of gas as an external loading. 

One should treat the integrand as a function of independent variables xu ^2, xi\ xi, t taking into 
account ( 12.3 ) and ( 12.4 ). 

12.4. Equation derivation principle 

To derive necessary equations, we use a very general technique known as Pontryaghin principle of 
maximum [ 18 ], [ 16 ]. This technique was developed to solve problems of optimal control, and 
both Hamiltonian and Lagrangian mechanics can be treated as special cases of it. 

The goal of solving a problem of optimal control is to find n state variables x\it\.,^n{t) and r control 
variables mi(0,..,W/<0 within a specific interval of variation of the argument fo ^ ^ ^ to minimize a 
given functional: 
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(12.8) ^=J/o(^,-. 

on the stipulation that the state equations be satisfied: 

( 12.9) ^' = /(^v,^«,«i,-,^r.O 0*=lv.,n). 
The stroke denotes the differentiation with respect to L 

Let's introduce n conjugate variables N\(t\..,Nn(t) and build die Hamiltonian function: 

(12.10) H(.v.^.,N,,..,N„,/.,,..,/.,,/)=Vo+Z^^^^ 

1=1 

The Pontryaghin principle of maximum states: the optimal control that minimizes the functional ( 
12.8) is formed by controls wi,..,Wr that maximize the Hamiltonian function for each t from the 
interval fo ^ ^1. 

Also, the following differential relationships must hold: 

(12.11) n; = ~ (i=l,..,n); 

(12.12) ^'=-|^ (/=l,..,n). 

Relationships ( 12.11 ),( 2.18 ) are canonical Euler equations for the problem being solved. The 
system of equations ( 2.18 ) is identical with the system of state equations ( 12.9 ). 

12.5. Formulation of the problem as an optimal control problem 

The problem of minimization of the functional ( 12.7 ) can be treated as an optimal control problem 

if we treat xi, X2 as state variables and xi xi' as control variables. Variables Xi \ X2 are assumed to be 
independent from jci, X2. Let's denote them as wi, wa- 
In this case the state equations are the following 
( 12.13) xx^uu 

( 12.14) X2'=M2. 

Thus having introduced the notation of the Pontryaghin principle, we have: 

( 12.15 ) /o= (ef +2v£,e, +e2)-PX(a- 0^1^2(1+ £2); 

(12.16) /,=«,; 

(12.17) /2 = «2; 
Notice that 

(12.18) e,=5'-l; 



(12.19) 



1 +"2 



Introducing conjugate variables M and A/2, we build the Hamiltonian function: 

( 1 2.20 ) H= (ef + 2ve,e2 + )+ ptifi - f)x^u^ (1+ ) + iV.u, + iV^u^ . 

2(1 V ) 
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For our further computations we need derivatives of the Hamiltonian function with respect to all 
variables. 

Derivatives with respect to the state variables: 

(12.21) =p^(^a-t)u,(l + E,); 

(12.22) ■T^ = 0- 

Derivatives with respect to the control variables: 



dH Eh 



u 



( 12.23 ) — = J (e, + ve^) ^ + ; 

1-V s 

( 12.24 ) — - = ^(e, +ve^)^+N, + px(a-t)x,il+E,). 

OUj 1 — V s 



Derivatives with respect to the conjugate variables: 

dH 

f2 



(12.25) TTr = u,; 



(12.26) — = 



According to the maximum principle of Pontryaghin, the optimal control that minimizes the 
functional ( 12.15 ) is achieved by the controls uu U2 that maximize the Hamiltonian function for 
each t from the interval 0<t< a-\-b. There are no limitations (in our particular case), so at least we 
should meet the conditions of stationarity of the Hamiltonian function ( 12.20 ) widi respect to uu 
M2. These conditions can be obtained by equating ( 12.23 ) and ( 12.24 ) to zero! 

(12.27) _^(e,+ve,)i^l: = A^,; 

1-v J 

( 12.28 ) ^(e, +ve^)!^^ = N,+pxia-t)x,il + E,) 

1 — V s 

These conditions are algebraic equations that establish relationships between the functions of~ 
interest. 

Also, the system of differential equations ( 12.11 ), ( 2.18 ) should be fulfilled: 

dN 

(12.29) ^ = .p^^a-t)u,a + E,); 

at 

(12.30) ^ = 0; 

dt 

(12.31) 4^ = «n 

dt 

(12.32) ^ = M,. 

dt ^ 
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Besides the equations ( 12.27 ) to ( 1232 ), four boundary conditions should hold according to the 
number of differential equations of first order in the system: 

( 12.33 ) JCi(0)=0; 

(12.34) JC2(0)=0. 

These conditions are obvious from Fig. 12.2. The same figure gives us: 

(12.35) ^|,,=-,(0)=0. 

at 

If we remember that jci(O) = 0, then the equation ( 12.24 ) gives us the following: 

(12.36) A^2(0) = 0. 

The fourth condition is formulated for the point t = a^-b. 
( 12.37) JCi(a+Z?) = a+/?-f A; 

where A is the forced elongation of the elementary cell. 

If we solve a problem with a given linear elongation of the airbag, then 

(12.38) G,(a + fc)=^; 

h 

so, taking into account the fact that in the vicinity of r = a+b the equahty u\=^s' holds, the equation ( 
12.27 ) gives: 

(12.39) N,ia + b) = q. 

12.6. Closed analytical solution of the problem 

It is intuitively obvious that in the interval 0 < r < a the integral curve is a circle, for example: 

(12.40) Xi=/?sin— ; 

R 



(12.41) Xj=/? 



l-cos — 



where R and k are arbitrary constants. 

This curve fulfils the initial conditions ( 12.33 ), ( 12.34 ). 

The solution for the conjugate variable is also obvious: 

(12.42) A^2 = 0. 

It satisfies both the equation ( 12.30 ) and the initial condition ( 12.36 ). 

In order to meet the requirement of satisfying the equations ( 12.31 ), ( 12.32 ), we derive the 
following expressions for the functions mi, M2 at the interval 0<t<a: 

(12.43) M, =jkcos— ; 

R 

(12.44) =ifesin— . 

Now we have from the formula ( 12.19 ): 

(12.45) s' = k; 
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and then from the formula ( 12.18 ): 

(12.46) ei=it-l. 

From the formula ( 12.27 ) we derive: 

Eh ^ . kt 
— -l + vejcos- 
l-v'^ R 



(12.47) N, =-^(jk-l + v£jcos^. 



By differentiating the expression of M and substituting the result to ( 12.29), we have the following 
after some simple transformations: 

Fh 

( 1 2.48 ) — - 1 + ve, ) = pR{\ + e J . 

The physical meaning of this formula becomes clear if we transform it as follows using the formulae 
( 12.5 ) and ( 12.46 ): 

(12.49) 0.=^^^^. 

' h 

From the formula ( 12.48 ) we can determine the arbitrary constant k: 

(12.50) ,=M±ii)M)^l_ 

Eh ^ 

The formula ( 12.47 ) becomes as follows: 

(12.51) A^, =p/?(l + e2)cos— . 

R 

Thus we have found a partial solution of the system of differential and algebraic equations ( 12.27 ) 
to ( 12.32 ) that satisfies the initial conditions ( 12.33 ), ( 12.34 ), ( 12.36 ) and contains one arbitrary 
constant R. 

The value of this arbitrary constant should be chosen using ( 12.37 ) or ( 12.39 ). Though, to do that, 
we need to extend the solution of the system of equations to the interval a<t^ a+b. 
Because at a < r < a-^b 

(12.52) X(«-0 = 0; 

the equations ( 12.28 ) and ( 12.29 ) at this interval have the following form: 

(12.53) _^(e +ve,)% = A^,; 
1-v s 

(12.54) ^ = 0. 

dt 

The following functions satisfy the equations ( 12.27 ), ( 12.53 ), ( 12.54 ), ( 12.30 ), ( 12.31 ), ( 
12.32): 

(12.55) j:, =/?sina + itj,(r-a); 

(12.56) =/?(l-cosa); 

( 12.57 ) A^, = p/?(l + £2 )cos a ; 

(12.58) N2 = 0\ 
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where 

(12.59) a = — ; 

R 

(12.60) K=^—^ — ^cosa + l-ve,. 

"Eh ' 

The formulae ( 12.28 ), ( 12.29 ) yield 

(12.61) ui = kt,; "2 = 0; s' = kb;Ei = kt, - \; 

and these expressions enable us to check the dependencies ( 12.55 ), ( 12.58 ) by the direct 
substitution to the respective equations. Also, the integral curve maintains its continuity in the 
vicinity of the point t = a. 

The stresses in the airbag material can be easily found from the formulae ( 12.5 ),( 12.6 ): 
(.2.62) ot=^^^: 



(12.63) a;=-^^&±^v + &.; 



(12.64) 



' h 

(12.65) cl ^^^^j^vcosa-^EE,; 

where is the stress along the guide in the non-glued area; 

al is the stress along the generatrix in the non-glued area; 
af is the stress along the guide in the glued area; 
02 is the stress along the generatrix in the glued area. 

The force that tears the sheets apart on die gluing boundary can be found from the conditions of 
equilibrium: 

( 12.66) r = /7/?(l-he2)sina. 

The aAitrary constant R can be found from the equality ( 12.39 ) in the following form: 

(12.67) = p/?(l + e2)cosa. 

To calculate the elongation A, we can use the formulae ( 12.37 ) and ( 12.55 ): 

(12.68) S = /gsina+ ^^^^'^^^^^^"^'K osa-a-v£2fc. 

Eh ^ 

If the value of 8 is specified, then the arbitrary constant R should be determined from the equation ( 
12.68 ) and the linearly distributed tension force from the formula ( 12.67 ). 

All computations above are made with the assumption that £2 is a given value. If the airbag is not 
prevented from deforming along its guide, then this value can be found from the condition of 
equilibrium of forces applied to the cross-section of the airbag: 

(12.69) Sp = olah + G'^bh; 
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where S is the area of the cross-section of the hollow filled with gas included in one elementary cell. 
By expanding the formula ( 12.69 ), we obtain: 

( 1 2.70 ) -^-^ (a - cos a sin a) = pR{a + b cos aXl + £2 V + Eh{^ ^b)^^. 

Thus, one needs to solve the system of two nonlinear equations ( 12.67 ) and ( 12.70 ), or ( 12.68 ) 
and ( 12.70 ), with respect to /? h £2 in order to find these values. 

12.7. Approaching the reality 

In the preceding sections we have obtained an accurate solution of the problem that describes an 
infinitely long model No.l. In this model, the cylindrical hollows are perpendicular to the tension 
force appUed to the airbag. 

The real system is finite, and it contains a finite number of hollows filled with gas. Also, it has 
additional elastic parts and elastic elements that affect its behavior as a whole. 
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Fig. 12.3. Design model of the airbag 

Fig. 12.3 shows a design model of the airbag. This model includes a working area that contains n 
hollows with the effective width c, and an additional elastic part of the compliance z. 

Let the tension force for this airbag be Q, and the elongation be A. On the basis of data obtained 
from the investigation of the elementary cell, we can determine these values by the following 
formulae: 

(12.71) Q = 2qc; 

(12.72) A = 25n + ez. 

12.8. The set of computational formulae 

Source data: 

h is the film thickness; 

E is the elasticity modulus of the film; 

V is the Poisson ratio of the film; 

a is the half- width of the non-glued band (Fig. 12.1); 

b is the half-width of the glued band (Fig. 12.1); 

c is the effective width of the bag (Fig. 12.3); 

n is the number of inflatable hollows; 
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p is the pressure within the bag; 

A is the elongation of the airbag comparing to its non-deformed length under the zero 
intemal gauge pressure. 

To analyze the airbag, one needs to compose two nonlinear algebraic equations: 



( 12.73) 

(12.74) 
where 
( 12.75 ) 



In 



/?sina + 



p/?fc(l + £,)(l-V^) 

Eh 



cosa-a-vfce. 



+ 2zcpR(l + 82 ) cos a = A ; 



«n2 

^^y- (a -cos a sin a)- pRv{\ ^E^X^ + b cos a) -Eh{a^b)E2 =0; 



^^ pa(l + £j(l-v^) ^ {l-VE,)a 
Eh R 



These equations should be solved (numerically) with respect to R and £2 and thus find the values. 
Next, using the formulae ( 12.62 ) - ( 12.66 ), ( 12.76 ), ( 12.77 ), one can find all values of interest. 

(12.76) d = 2R{l-cosa); 

( 12.77 ) Q = 2cpR{l + £2 )cos a ; 

where d is the thickness of the inflated hollow (Fig. 12.3); 
Q is the tension force. 

12.9. Behavior of a separate hollow in a stiff airbag 

Peculiarities of the behavior of the structure can be revealed by analyzing a separate hollow 
contained in an aiibag made of an infinitely stiff material (Fig. 12.1). 




Fig. 12.4. The model of a separate hollow 

Let's assume: n = I; b = 0; E = 00; c = I. 
Then the principal formulae will look like: 



( 12.78 ) 


a 


( 12.79 ) 


= 2/?(l-cosa) 


( 12.80 ) 


/ = 2/fsina; 


( 12.81 ) 


A = 2(/?sina-a 


( 12.82 ) 


= 2p/?cosa ; 
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The latter foUows from ( 12.69 ). 

These formulae contain the only arbitrary constant R that can be determined either from the equation 
12.81 ) or from the equation ( 12.82 ) depending on the source data. 

Let's introduce dimensionless values: 



12.83 ) 
12.84 ) 

12.85 ) 



r = pRsina; 



o_pR 
h 



pR 
2h 



''j Sin 2a ^ 



V 



2a 



12.86 ) 



1i 



a = — ;/= — ; A = — ; q = — —; r = 
2a la la Ipa Ipa 



Ipa 



^2 



2 pa 



Then the formulae ( 12.79 ) - ( 12.85 ) can be written in the form: 
12.87 ) 



T 1-cosa 
a = 



12.88 ) 
12.89 ) 
12.90) 
12.91 ) 
12.92) 
12.93) 



/ = 



a 
sin a 



a 



-r Sin a 

A = 

a 

_ cos a 



__ sma 
2a 

1 

' 4al 



-1; 



sin 2a ^ 
2a 



Here the arbitrary constant is the angle a. 

The dependencies are shown graphically on Fig. 12.5 - Fig. 12.7. In all dependencies the argument 
is the dimensionless elongation A . The limits of variation of this elongation, provided the force q 
is positive (extending the airbag), can be found from the inequalities: 



(12.94) 



--1<A<0. 

K 
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-0.4 -035 -0.3 -0.25 -0.2 -0.15 -0.1 -0.05 0 
Dimensionless elongation - delta 

dimensionless force q 

dimensionless force r 



2/pi - 1 

Fig. 12.5. Dependence of dimensionless forces ^and r on the dimensionless elongation A 




-0.4 -035 -03 -0.25 -0.2 -0.15 "O.l -0.05 0 
Dimensionless elongation - delta 

dimensionless tlength 1 

dimensionless thickness d 



2/pi - 1 

Fig. 12.6. Dependence of relative sizes of the hollow / and d 
on the dimensionless elongation A 
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-0.4 -0.35 -0.3 -0.25 -0.2 -0.15 -0.1 -0.05 
Dimensionless elongation - delta 
dimensionless stress sigmal 

dimensionless stress sigma2 

2/pi-l 



Fig. 12.7. Dependence of dimensionless stresses and on the dimensionless elongation A 
12.10. Behavior of the real airbag 

Let's take a look at a real airbag for the "Satum" car (see Fig. 1 1. land [ 21 ]). This airbag has four 
inclined cylindrical pockets in the driver area and another pocket of the same kind in the passenger 
area. The hollows are about 24 cm long and 10.9 cm wide. Glued stripes between them are about 1 
cm wide. The material is a synthetic fabric with impregnation. 

We assume the following source data: actual values fit = 5.45 cm; b = 0.5 cm; n = 5; intuitive values 
c=l2cm;h = 0.01 cm; E = 4000 kg^cm^; v = 0.39; z = 0.25 cm/kgf. Using the results of Section 
12.8, let's build curves of deformation of the side impact airbag. The curves are shown on Fig. 12.8. 
The comparison of these data with experimental results on Fig. 12.9 shows that the general behavior 
of the dependencies and the order of magnitude of values is described correctly by the approximate 
theory. 
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0.4 0.5 0.6 
pressure kg&sq.cm 



delta = 0 
delta = -2 
delta = -4 
delta = -6 
delta = -8 
delta = -10 
delta = -12 
delta = -14 



Fig. 12.8. Dependence of the longitudinal tension force on the internal pressure and elongation 




pressure kg&sqxm 

d=-14cm 

d=-12cm 

d=-10cm 

d=-8cin 

d=-6cin 

d=-4 cm 

d=-2 cm 



Fig. 12.9. Experimental data from the paper [ 21 ] 
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Theoretical results have been compared with experimental data. In spite of the model being very 
different from the real structure, results of the comparison can be treated as quite satisfactory. 

13. AIRBAG WITH HOLLOWS AT AN ANGLE TO THE TENSION 
13.L Introduction 

In this chapter we will investigation a simplified model of a side impact protector airbag that keep 
the essential behavior of the real object of interest. 

The chapter 12 deals with a model in which there are cylindrical hollows filled with gas and 
perpendicular to the longitudinal axis of the airbag. The present chapter deals with a model in which 
the hollows are at some angle to the longitudinal direction of the airbag. As in the previous case, the 
goal of the study is to find out dependencies between die gas pressure, the airbag tension and its 
elongation. 

The chapter derives differential equations, gives solutions to them in a closed form, and builds sets 
of curves for a non-deformable film and for the airbag model with which an experiment has been 
run. 

13.2. Design model 

As in the chapter 12, we assume the system of interest to include an infinite sequence of cylindrical 
hollows spaced evenly, each hollow being infinitely long. The system is built out of two infinite film 
pieces glued with stripes 2b wide. Between each couple of adyacent glued stripes there is a non- 
glued gap 2a wide that makes a cylindrical hollow filled with gas under the pressure p. Thin film 
that the bag is made of has the thickness h, elasticity modulus E, Poisson ratio v. We assume the 
film not to be resistant to flexion. 




Fig. 13.1. A cut from a system having a sequence of infinitely long cyhndrical hollows 
perpendicular to the lengthwise direction of the airbag 

Fig. 12.1 shows a cut fragment of an infinite system with the length 1 along the generatrix. The cut 
includes two cyhndrical hollows and their adjacent stripes of gluing of the width b. The sizes a, b 
and the unit size refer to the non-deformed state of the film. In the course of loading, the fihn 
deforms and the sizes change. 

As one can see from the above said, the shape of the bag is drawn in space by a rectilinear generatrix 
moving along a curviUnear guide. 

The cloth of the bag is tensioned at an angle to the guides, the component of the linear load 
perpendicular to the guide is 2q, and the component acting along the generatrix is 2w. 

This system has certain properties of synmietry, and we can analyze only a repeatable elementary 
cell instead of the whole system [ 19 ]. 
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Fig. 13.2. An elementary cell of the system from Fig. 12.1 

Fig. 12.2 shows fihn fragments belonging to the elementary cell with bold lines in cross-section. 
One should note that the elementary cell includes only one of the glued films and only one quarter of 
the gas that fills one hollow. The elementary cell accepts a half of the linear load that tensions the 
airbag. 

133. Mathematical formulation 

Let's refer the elementary cell to a Cartesian coordinate system jci0jc2 as shown on Fig. 12.2. We 
introduce into consideration the parameter t which has the following physical meaning: it is the arc 
distance along the guide line measured on the non-deformed film surface. The reference point of the 
measurement is the middle of the non-glued area of the film (Fig. 12.2). 
The parameter may change within the following limits: 
(13.1) 0<r<a+fr; 

where a is the half-width of the non-glued airbag area in the non-deformed state; 

b is the half-width of the glued airbag area in the non-deformed state. 

The value t = a corresponds to the boundary between the glued and non-glued areas. 

As the airbag is a cylindrical body with a rectilinear generatrix, its deformed state is completely 
defined by the section of the guide line described above. 

Let this line be specified parametrically using coordinates x\(t\ JC2(0, x^it) that depend on the 
parameter t. 

The displacement method will be used to solve the problem. During the derivation of formulae we 
will introduce new functions, all depending on the parameter t. 

As all variables are fiinctions of the only parameter f, its notation as an argument of the functions 
will be omitted. The differentiation with respect to this argument will be denoted by a stroke. 
Apparently, the elementary arc length in the deformed state is described by the formula: 
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( 13.2 ) ds = dtylixlf+(x^y+(x;f . 



The relative elongation of the film along the guide curve Ei can be described by the following 
formula: 

( 13.3) £1 = — - — = 5 -1; 

at 

where 

(13.4) s^^{x[f+{x\f^{x[f . 

The relative elongation along the generatrix will be denoted as £2. Note that this is a constant that 
does not depend on the parameter t. 

The distortion Y in this case is described by the following formula [ 17 ]: 

(13.5) Y = 4. 

s 

To define the elastic potential energy of the system and the stresses in the film, we need physical 
equations: 

( 13.6) Gi r(e,+V£2); 

1-v 

(13.7) 02=— ^(vEi+eJ; 

1-v 

(13.8) x= ^ 



2(l + v)'' 

where d and <52 are normal stresses along the guide and the generatrix respectively.; 
X is the tangential stress. 

13.4. Full energy of the system 

Here we use a variation approach to the displacement method solution. This approach gives us a 
good formalization of this solution phase, helps choose a consistent system of variables, formulate 
natural boundary conditions, pose the dual problem if necessary. Also, there is no need to compose 
equilibrium equations looking at figures. 

Let's choose a system of fiinctions and constants to completely define the deformed state of the 
system. These fiinctions can be xi, JC2, and the constant £2. 

The full energy of the system W, that includes the energy of loading also, is now described by the 
following fimctional: 



a+b 

( 13.9 ) W = \ 



Eh ( "> ^ o\ Eh <■) / 

[e\ + 2V£j£2 + £^ j+ -r r f " PX(fl " t)x,X^ (1 + E^ ) 



2(1-V^)' ' ' ' 4(1 + v) 



dt\ 



where x(z) is the Heaviside fiinction equal to 1 at z > 0 and 0 at z < 0. 

The first term in the brackets is the elastic energy of the airbag material, and the second term is the 
energy of gas as an extemal loading. 

One should treat the integrand as a function of independent variables xu X2, x\\ X2, xz\ t taking into 
account ( 12.3 ) and ( 12.4 ). 
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13.5. Equation derivation principle 

To derive necessary equations, we use a very general technique known as Pontryagin principle of 
maximum [ 18 ], [ 16 ]. This technique was developed to solve problems of optimal control, and 
both Hamiltonian and Lagrangian mechanics can be treated as special cases of it. 

The goal of solving a problem of optimal control is to find n state variables xi(t),,.^n(t) and r control 
variables u\(t),..,u/^t) within a specific interval of variation of the argument tQ< t <ti to minimize a 
given functional: 

( 13.10 ) >^ = j/o(>ir., ^.,^1,.-, ^ry^d^ ; 

on the stipulation that the state equations be satisfied: 

( 13.11 ) ^' = /(>iv.,^„,^i,..,/^r,0 0*=l»-.»'i). 
The stroke denotes the differentiation with respect to t. 

Let's introduce n conjugate variables M(0» -^nCO and build the Hamiltonian function: 

(13.12) H(:s..^.,N,,..,N„,//,,..,/.,,/)=-/o+E^^y:- 

1=1 

The Pontryaghin principle of maximum states: the optimal control that minimizes the functional ( 
12.8) is formed by controls Mi,..,Mr that maximize the Hamiltonian function for each t from the 
interval ro^ ^ ^ ^i- 

Also, the following differential relationships must hold: 

(13.13) H = -T- 0=1.-.«)J 

(13.14) x;=^ (/=l,..,n). 

Relationships ( 12.11 ),( 2.18 ) are canonical Euler equations for the problem being solved. The 
system of equations ( 2.18 ) is identical with the system of state equations ( 12.9 ). 

13.6. Formulation of the problem as an optimal control problem 

The problem of minimization of the functional ( 12.7 ) can be treated as an optimal control problem 
if we treat xu x^, as state variables and jci', x^, jca' as control variables. Variables jci', X2^ jcs' are 
assumed to be independent from x\, jc2, X3. Let's denote them as wi, M2, W3. 
In this case the state equations are the following 

(13.15) xr=wi; 

(13.16) JC2'=M2; 

( 13.17) X3'=W3. 

Thus having introduced the notation of the Pontryagin principle, we have: 

( 13.18 ) /o = ^^^^^^,^ (£f +2v£,£, +E2)+^j^y' -PX(«-0^iW2(l + e2); 

(13.19) /i = wi; 

(13.20) /2 = M2; 
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(13.21) /3 = W3. 

One should note that 

(13.22) ei=5'-l; 

(13.23) y = %. 

s 

(13.24) s'^yjuf+ul+ul ; 

Introducing conjugate variables M, A'iz and N^,, we build the Hamiltonian function: 
( 13.25 ) 

H = " 2(i^v2^ fe + 2veie2 + )+ y' + /7X(a - t)x,u, (1 + ) + A^iW^ + iV2"2 + A^3«3 • 

For our further computations we need derivatives of the Hamiltonian function with respect to all 
variables. 

Derivatives with respect to the state variables: 

(13.26) =p5^(a-OM3(l + e,); 

(13.27) -?^ = 0; 



dx 



2 



(13.28) ^ = 0. 

Derivatives with respect to the control variables: 

( 1 3.29 ) — = - ^— J (e, + ve^ ) -L + + ; 

^" ' " s 2(1 + vX^) 

(13.30) ^ = .-3(e,+ve,)^^ + ^g^ + /V,+ra(a-.)..,aH-e,); 

(13.31) ^ = .^(e,..ve.)%+^*??!i^+N3-. ^^l- 





1 








Eh 




1 




a// 




Eh 


auj 


1 





5' 2(1 + vX/)' ' 2(1 + v>' 



Derivatives with respect to the conjugate variables: 

(13.32) §--u,, 

(13.33) ^ = "2; 
(.3.34) ^ = „3. 

According to the maximum principle of Pontryaghin, the optimal control that minimizes the 
functional ( 12.15 ) is achieved by the controls mi, W2, U3 that maximize the Hamiltonian function for 
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each t from the interval 0< t< a+b. There are no limitations in our particular case, so at least we 
should meet the conditions of stationarity of the Hamiltonian function ( 12.20 ) with respect to ui, 
U2, M3.. These conditions can be obtained by equating ( 12.23 ), ( 12.24 ) and ( 13.31 ) to zero: 



( 13.35) 



(13.36) 



( 13.37 ) 



Eh 


fs' 


-l + VEj 


ul ] 




1 + v 




1-v 




s 


Eh 


f / 


-l + VEj 


2 ^ 


ill 


1 + v 




1-V 


2(/)^J 


s 


Eh 


(s' 


-l+ve2 


2 > 




1 + v 


\ 


1-v 




f 

s 



Ehu. 



2(1 + vX/)^ 



These conditions are algebraic equations that establish relationships between the fimctions of 
interest. 

Also, the system of differential equations ( 12.11 ),( 2.18 ) should be fulfilled: 



( 13.38 ) 
( 13.39 ) 
(13.40) 
( 13.41 ) 
( 13.42 ) 
( 13.43 ) 



dt 

dt 

dN. 



= -p%{.a-t)u^(\ + z^)\ 



= 0 



dt 



3 _ 



= 0; 



dt 



= u 



1 » 



2 _ 



dt 



2 ' 



dx 



dt 



Ideally, we should have determined uu W2, M3 from the equations ( 12.27 ) - ( 13.37 ) and substituted 
those to the differential equations ( 12.27 ) - ( 12.32 ). Then these equations would contain only six 
unknown functions. 

Besides the equations ( 12.27 ) - ( 12.32 ), we should meet six boundary conditions according to the 
number of equations of first order in the system. 

Conditions at the beginning of the integration interval are obvious from Fig. 12.2: 

(13.44) xi(0)=0; 

( 13.45 ) JC2(0)=0; 

(13.46) X3(0)=0. 

Conditions at the end of the integration interval are defined by the requirement that the elongation of 
the elementary cell in a skewed direction (at the angle (p as can be seen on Fig. 13.3) be equal to a 
given value of 5. The angle (p is also known. It can be measured on the airbag directly. 
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Fig. 13.3. A skewed cut from the system having a sequence of infinitely long cylindrical hollows 
The manifold of final states [ 16 ] is described by the following fomiula in this case: 
( 13.47) G = 4xl+xl -(a + fcyi + tan'9"-5 = 0; 

where 

( 13.48 ) =(a + fcXl + £2)ta^^9 + ^3- 

Three more boundary conditions are the trans versality conditions [ 16 ]: 



( 13.49 ) 



CtC; 



They look as follows: 



(13.50) 

( 13.51 ) 
(13.52) 



___£fL_- 



= 0 



(r = a+fo; 1= 1,2,3). 



{t - a+b); 



N2 = 0 



9^4 



I 2 2 



= 0 



(t = a-^b); 



where ^ is a new unknown constant. To determine the latter, we add the equation ( 13.47 ) to the 
system ( 13.50 ) - ( 13.52 ), in the following form: 



( 13.53 ) 



yjx^ +JC4 -(a + feyi + tan^(p-"A = 0 



(t = a-^-b). 



13.7. Closed solution of the problem 

The analytical closed solution has been guessed. It is given in the table that follows, with 
undetermined constants /?, fc, ^, ^1, T. 





0<f<a 


a<t<b 


( 13.54 
) 


X. =/fsin — 
' R 




jCj =/?sm — + fci(r-a) 
R 


( 13.55 
) 


Xj = R\ 1-cos — 




=R\ 1-cos — 


( 13.56 
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) 

( 13.57 
) 

^ 1 

^ Ij.JO 

) 

( 13.59 
) 


kt 

Ni = p/?cos— (l + Ej) 
R 

=0 


= p/? cos (1+82) 
^2 =0 
=7 


The auxiliary functions follow: 




0<t<a 




( 13.60 
( 13.61 
( 13.62 
( 13.63 
( 13.64 

( 13.65 

( 13.66 
( 13.67 

( 13.68 


kt 

u, =A:cos — 
' R 

, . kt 
' R 

"3 =C 

x= 


M, =/:, 
M2 = 0 

"3 = Ci 

ei=7fc'+Ci -1 
ai-j_^2(A'+Ci l + ve2) 


2(l+v)Vit' + C' 


2(l + v)V^,^+C, 



A direct substitution makes it sure that the differential equations ( 12.29 ) - ( 12.32 ) hold. Also, the 
initial conditions ( 13.44 ) - ( 12.34 ) and the condition ( 13.51 ) hold. Besides, the integral curves 
are continuous in the vicinity of the point t = a. 



To comply with the constraint equations ( 12.27 ) - ( 13.37 ) and boundary conditions ( 13.50 ), ( 
13.52 ), ( 13.53 ), we should choose the constants /?, k, ku C> Ci» q correctly. 
A system of algebraic equations to determine these constants follows: 

( 13.69) p/t , ^ =/?/?(l + £j; 
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( 13.70) 
( 13.71 ) 
( 13.72 ) 

( 13.73 ) 

( 13.74) 
( 13.75 ) 



ph 



= /?/?(l + £2)cosa; 

Eh^ 



= T- 



= T- 



2(l + v)(^'+C')' 



.[j^^ 2(l + v)(^,^+Cf)' 



p/?(l + )cos a = (/? sin a + ^,&) ; 
= [(a + *)(1 + e 2 ) tan 9 + aC + ] ; 



( 13.76 ) 

( 13.77 ) 

( 13.78 ) 
( 13.79) 



Z) = (a+fr)7l + tan^ + 5. 

)een introduced t< 



Here the following notation has been introduced to reduce the length of the expressions: 

f <-- ^ 



P = 



1+v 



1-v 



pi = 



l + v 



l-v 



£> = ^[(a + fc)(l + e2)tan9+a^ + i^if +(/?sina + ikife)^ ; 



a = — . 



The equations ( 13.69 ) and ( 13.71 ) are the equations ( 12.27 ) and ( 13.37 ) written for the interval 
0 < r < a, and the equations ( 13.70 ) and ( 13.72 ) are the equations ( 12.27 ) and ( 13.37 ) written 
for the interval a<t<b. The equations ( 13.73 ), ( 13.74 ), ( 13.75 ) are equivalent to the equations ( 
13.50 ), ( 13.52 ), ( 13.53 ). The equation ( 12.28 ) is not used here because it is an identity as long 
as the equation ( 12.27 ) holds. 

So, we have seven equations to determine seven constants. 

All calculations above use the assumption that £2 is a given value. If the airbag is not constrained 
against its deformation along the generatrix, this value can be found from an equation of equilibrium 
of forces applied to the airbag cross-section: 



( 13.80) 



Sp = (5lah + <5\bh; 



where S is the area of the cross-section of the hollow filled with gas per one elementary cell; 
^2 = ^zlo^sa ' <^2 = <^2\a<,ia*b ^ according to the formula ( 13.67 ). 
We can derive the following by expanding the formula ( 12.69 ): 



( 13.81) 



n2 

(a - cos a sin a) = pR{a + b cos aXl -»- £2 )v -h Eh{a + b)E..^. 
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Thus, in order to determine 62, we should add the equation ( 12.70 ) to the system ( 13.69 ) - ( 
13.75 ) and solve all these equations as a single system. 

For the computation purpose, it is reasonable to use another system equivalent to the equations ( 
13.69) and (13.71): 



( 13.82) 



ph = Jp'R'{l + ^J + 



T - 



Ehtj 

2(l + vXik' +^') 



( 13.83 ) 

use an equivalent system instead of the equations ( 13.70 ) and ( 13.72 ): 
( 13.84 ) 



p,/i = 



p^R^il + E^y cos^a + 



T- 



2(l + vX/:,^+Cf). 



|pMl.eJcosa = r-^^^^^^; 



( 13.85 ) 

also, use an equivalent system instead of the equations ( 13.73 ) and ( 13.74 ): 



( 13.86) 
( 13.87 ) 



^ = ^p2/j2(i+eJ2cos'a+r' ; 
(a + fe)(l + e2)tan(p+aC+feCi = 



T[Rs.\na-\-kib) 
/?/?(l + e2)cosa 



In this system, an unknown value can be expressed via other values using the direct formula ( 
13.86 ). Thus the system is reduced by one equation. 

13.8. Finding unknown constants by a direct method 

In order to find the unknown constants R, k, k\, C„ £2, we can find the point of minimum of the 
function ( 12.7 ) in the space of these unknown constants observing the limitation equality (13.75 ), 
instead of solving the system of nonlinear equations ( 13.69 ) - ( 13.75 ), ( 12.70 ). 

To do that, we substitute the expressions of x\, xi, Z\, y from the formulae ( 13.54 ), ( 13.55 ), ( 
13.64 ), ( 13.65 ) to ( 12.7 ) and then take the define integral thus representing Was a function of the 
desired parameters. A standard procedure of Mathcad-2000 determines die needed constants 
numerically. Let's do this exercise: 

( 13.88) 

(,le+C-i!M'Jk'+C,'-l^,+El ^ C 

1-v 2(/fe'+C') 



W=- 



Eha 



2(1 + v) 
Ehb 



2(1 + v) 



(Vfc.^+;?-lf+2v(Vfc.^+C?-l)k,+E^ ^ C 



1-v 



pkRaf, R . 2ka\ 

1 sm 

2 2fa2 R J 



Note that the accuracy of determining this fiinction's minimum is not always satisfactory. Therefore 
the further investigation is based on the solution of the system of nonlinear algebraic equations. 
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13.9. Changing to the real system 

In the previous sections we have obtained an exact solution of the problem regarding a simplified 
infinitely long model No.2. In this model its cylindrical hollows are placed at an angle to the tension 
force that stretches the airbag. 





............. 


1 

y. 


// // /ljp_ 


11 1 J ^ 


m 



A- A 



OC^XIXD ^AAA AAr^ 



Fig. 13.4. A design model of the airbag 

The real system has finite sizes and a finite number of hollows filled with gas. Also, it has additional 
elastic parts and elements that affect its behavior as a whole. 

Fig. 12.3 shows a design model of the airbag. This diagram includes a working area with a certain 
number of hollows of the effective width c, and an additional elastic part with the compHance z. The 
inclination of the hollows is described by the angle (p as shown on the figure. 

Let the tension force be Q and elongation A for this airbag. Referring to data obtained for the 
elementary cell model, we can determine these values by the formulae: 

( 13.89 ) Q = 2^cVl + tan'(p ; 

(13.90) A = 25n + ez. 
13.10. Procedure of analysis 

It is reconmiended to obtain all necessary dependencies for a real airbag as diagrams (graphs). The 
convenient parameter for the graphs is the elongation of the elementary cell 5. The lower limit of 
this value is negative: this occurs when the pressure in the hollow tends to zero and there is no 
tension force applied to the airbag. The lower elongation limit can be expressed by the formula: 

( 13.91 ) =(a-ffc)(V4^+tan>-Vl + tan>); 

^ 2a-hnb 

where q = —7 r . 

7c(a + fc) 

The upper limit can be 0 in most cases. 
The source data for the problem are: 

hj the thickness of the film; 

E, the elasticity modulus of the film; 

V, the Poisson ratio of the film; 

a, the half-width of the non-glued stripe (Fig. 12,1); 
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b, the half-width of the glued stripe (Fig. 12.1); 

c, the effective width of the airbag (Fig. 12.3); 
n, the number of the inflatable hollows; 

/?, the pressure in the airbag. 

For each value of the parameter 5 we should solve the compatible system of seven nonlinear 
equations ( 13.75 ), ( 12.70 ) - ( 13.85 ), ( 13.87 ) with respect to the variables /?, 82, ku ^, T. 
Then, we can find the stress in the material by the formulae ( 13.66 ) - ( 13.68 ), the tension force per 
unit length by the formulae (13.86 ), the full tension force in the airbag and its full elongation by the 
formulae ( 13.89 ) and ( 13.90 ). 

Also, we can determine the thickness of the inflated hollow (Fig. 12.3) 

(13.92) J = 2/?(l-cosa); 

and the force per unit length that tears the film sheets apart at their glued boundary: 

( 13.93) r = pR{l-hE^)s\na. 

13.11. Dependencies for the real airbag 

Here the real airbag will be a side impact head protector for the "Saturn" car [ 21 ]. This airbag has 
four inclined cylindrical hollows in the driver area and one more hollow in the passenger area. The 
hollows are about 25 cm long and 10.9 cm wide. The glued stripes between them are about 1 cm 
wide. The material is an impregnated synthetic fabric. 

The chapter 12 gives a set of curves built for the following source data: a = 5.45 cm; b = 0.5 cm; n = 
5; c = 12 cm; h = 0.02 cm; E = 4000 kg^cm^; v= 0.39; z = 0.2 cm/kgf. Fig. 12.8 dupUcates the data 
here for convenience of comparison. 

The fotogranmietry of the airbag (Fig. 11.1) shows that the four principal hollows are placed so that 
tan(p=0.528. Results of calculations done with the same source data and the inclination taken into 
account are shown on Fig. 13.6. 
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Fig. 13.5. Dependence of the lengthwise tension force in the airbag on the internal pressure and 

elongation, from chapter 12 

The comparison shows that the inclination of the hollows causes the lengthwise tension of the airbag 
to increase when the airbag shrinks Uttle, and vice versa: it causes the tension to reduce when the 
airbag' s length shrinks much. 
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Fig. 13.6. Dependence of the lengthwise tension force on the internal pressure, 
taking into account the inclination of the hollows 
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Fig. 13.7. Experimental data from [ 21 ] 

The model with the inclined hollows shows better correspondence to the reality, comparing to the 
model with the perpendicular hollows. 

Though, any comparison with experimental data is very relative. We don't know for sure the 
material thickness, nor its mechanical properties. The compliance factor of the system z is also 
unknown. It would be good to make a comparison between theoretical and experimental data after 
having measured the material thickness h, elasticity modulus £, Poisson ratio v, compliance factor z. 
As for the effective width of the aiibag c, we think it's acceptable to determine it by the formula: 
(13.94) c = /-2.4a; 

where / is the width of the area with the hollows (see Fig. 12.3); 

a is the half-width of the hollow (see Fig. 12.1 and Fig. 13.3). 

13.12. Effect of inclination of hollows 

Let us consider the dependence of the airbag tension force on the inclination angle of the hollows. 
To do that, we use the source data from Section 13.11, with the inclination angle tangent varied 
within the Umits 0 to 1.0. The analysis has been done for two levels of pressure: low 0.1 kg^cm^ and 
high l.Okg^cm^. 
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Fig. 13.8. Dependence of the aiibag tension force on the inclination of the hollows 

at the pressure 0.1 kgfi^cm^ 
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Fig. 13.9. Dependence of the airbag tension force on the inclination of the hollows 

at the pressure 1 kgf^cm^ 

Fig. 13.8 and Fig. 13.9 shows that the inclination of the hollows does not influence dramatically the 
mechanical behavior of the airbag. The principal tendency is that the airbag becomes more 
compliant in the lengthwise direction as the inclination angle increases. This tendency is not quite 
the case at a low gas pressure within the airbag. Thus, the inclination angle of the hollows can be 
chosen from engineering considerations, without referring to the airbag tension regulation. 

14. AXISYMMETRIC PROBLEM OF STRESS CONCENTRATION NEAR A ROUND 
GLUED AREA 

This chapter features a solution of an axisynmietrical problem for a simplest safety airbag made of 
two round sheets of film connected along their outline and also in a round area in their central parts. 
The following is to be investigated: the airbag' s shape under pressure, stress in the film including 
the round glued area, stress in the glued seam. Special attention is paid to the linear stress of the 
glued seam peeling and the peeling angle of the glued pieces of film. 

14.1. Design model 

The airbag under consideration consists of two round pieces of film of the radius b glued along their 
exterior contour and within a circle in their conunon central part of the radius a. The schematic 
model is shown on Fig. 14.1. This airbag has an annular hollow which is b - a wide in its initial 
non-deformed state. As the airbag is inflated, the width of this hollow decreases. 




Fig. 14.1. A sketch of the inflatable toric airbag. 
The glued area is grayed 

We assume the film to have the thickness of A, to be made of an isotropic material, and to have the 
elasticity modulus E and Poisson ratio v. 

The problem is to determine the stress and strain in the film when the hollow is filled with gas under 

pressure p above ambient. 

Let's assume the system may have large deformations and displacements. We assume also that folds 
are "smeared" (evenly distributed) over the material. 

Let the aiibag is referred to a cylindrical coordinate system r, 9, z so that the axis of revolution of 
the airbag be coincident with the z axis. As the airbag is a body of revolution, its shape is completely 
determined by a fragment of a curve that represents the meridional section of the shell. This curve is 
shown on Fig. 14.2 schematically. 
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Fig. 14.2. Cross-section of the toric airbag 

Let this line be specified parametrically using coordinates r(t\ z(t) depending on parameter t. 
Parameter t is the length of the line that connects the airbag' s center and the current point M when 
the airbag is not deformed. 

As the system is dissynmietric with respect to the plane z = 0, we will deal with only the lower semi- 
curve z < 0. Parameter t will then change within the following limits: 0<t<b, 

14.2. Physical relationships 

The following conventional physical relationships hold for the material in tension: 

(14.1) a, =-^(e,+V£j; 

1-v 

(14.2) <y2=Y^(vei+e2); 

El is an elastic relative elongation of the material in the meridional direction; 
£2 is an elastic relative elongation of the material in the circular direction; 
Ci is a meridional stress; 
O2 is a circular stress. 

14.3. Differential equations 

The displacement method is used to solve the problem. Its principal unknowns are r{t) h z{t). As all 
variables are functions of the only parameter r, its notation will be omitted further below. 

The problem posed is similar to that posed and solved in[ 20 ] in the part of derivation and solution 
of the appropriate differential equations. The difference is that there is no pressing die in this 
problem, and also boundary conditions are different. 

Let's give the system of differential equations omitting its derivation. 



Table 1. System of differential equations 



No folds 




Folds 


^,_iV,(A:-vA) 


U 




N 


r = — - — 
N 
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_,_N,{K-vA) 






N 




z = 

N 


^, _BA InrpN^iK -vA) 




^, _ InrpN^K 


'~ t N 




N 


n;=o 




Ar,=o 



Besides the principal unknown functions r and z, conjugate functions Nr and Nz feature in this 
system. The following notation has been introduced to simplify the writing: 



( 14.3 ) 


B=2itthE; 


(14.4) 




( 14.5 ) 




(14.6) 


B 


(14.7) 


A = — 1+ — . 



At A > 0 there are no folds, at A < 0 there are radial folds. 

In addition to the notation described in [ 20 ], the symbol ( 14.5 ) has been introduced and ro has 
been replaced by t. 

14.4. Boundary conditions 

Let's deal with the part of the airbag where the hollow is located: a<t<b. 

It is obvious that the z coordinate at the beginning and at the end of the integration interval will be 

zero. 

(14.8) z(a) = 0; 
( 14.9 ) z{b) = 0. 

As for the boundary conditions at the point t = a, note that at 0 < r < a the uniform tension of the 
central round area takes place, with the following equalities being true: 

( 14.10) £i = e2 = e; 

where e is an arbitrary constant. 
Thus 

(14.11) r = (l+e)f; 

(14.12) a,=c,=— . 

' ' 1-v 

Repeating the considerations from Section "Boundary conditions" of the paper [ 20 ], we come to a 
conclusion that the initial conditions for the unknown functions r and Nr can be expressed via the 
arbitrary constant e by the following formulae: 

(14.13) r(fl) = a(l + £); 
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( 14.14) 



InhEm 
1-v 



Also, 
(14.15) 



This is a natural condition for function Nr at the end of the interval. 

Function Nz is constant at the whole interval of integration. Its value is another arbitrary constant. 
Using tiie Runge - Kutta method for the integration of the system of differential equations, one 
should specify z{a) = 0 and then choose values of constants 8 and so that conditions ( 14.9 ) and ( 
14.15 ) hold at the end of the integration interval. 

14.5. Formulae for the desirable values 

In this chapter we are interested with the shape of the hollow in its inflated state and linear stresses 
of the glued seam peeling. 



The shape of the hollow can be described sufficiently well by three radii: large radius pi, medium 
radius p2, small radius p3, and the peeUng angle of the glued films a at the boundary of the central 
glued area (Fig. 14.3). Point M used here for the determination of all three radii corresponds to the 
ultimate value of the z coordinate. 

The radial coordinate of this point is found from the condition of the equality to zero of function A^o 
defined by formula ( 14.4 ): 



z 




1^ 



r 



Fig. 14.3. Parameters of the shape 



( 14.16) 




Thus 



(14.17) 



Pi 



ru -r{a) \ 




(14.19) p3=-imnz(0. 



Sine of angle a can be determined by the following formula: 
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(14.20) sin(a) = -^^. 

The linear peeling stress at the glued seam of the central area is determined by the formula: 

(14.21) T=_iW 

" 2iui 

The following formula determines the linear peeling stress of the exterior contour: 

(14.22) 7 -i^. 

The peeling angle at the contour is 90°. 
14.6. Dimensionless quantities 

Let us assume we are considering a curvilinear hollow that has the following half-width in its flat 
unstressed state (in plane): 

(14.23) 

2 

and its radius of rounding along the interior edge is a (cm. Fig. 14.2). 
Let's introduce a dimensionless parameter 

(14.24) 1] = -; 

a 

that describes the relative curvature of the hollow. 

Note that r| = 0 for a rectilinear hollow, and r| -> oo for a most curved hollow when a -> 0. 
Let's introduce dimensionless coefficients of shape: 

(14.25) e,=-BL; 

Po 

(14.26) e, =-^; 

Po 

/ 1 /I \ Q _ P3 . 



where po is the radius of the rectiUnear hollow for the case of an inextensible film, calculated by the 
formula: 

(14.28) Po= — . 

Let's introduce stress concentration coefficients: 

T 

( 14.29 ) -~ seam of the central round glued area, 

T 

( 14.30 ) K-~ for the seam of the exterior contour. 

^0 
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The reference level of the linear stress is assumed to be its value for the rectilinear hollow: 



( 14.31 ) 



_ 2pd 



n 



where p is an internal pressure above ambient in the airbag. 

This formula will be precise if the aiibag is made of an inextensible film. 
Similar derivations are made for the radial stress too. 
The stress concentration coefficients are determined in the following way: 

_a,(0) 



( 14.32 ) 

( 14.33 ) 

area, 

( 14.34 ) 



Co 

c. 
c.= 



"^0 



in interior points of the central round glued area, 
in the hollow region near the boundary of the central glued 

in the hollow region near the outer contour; 



where ai(0), (Ji(a) and 0\(b) are radial stresses in the film in interior points of the central round 
glued area, in the hollow region near the boundary of the central glued area, and in the hollow region 
near the outer contour, respectively. 

The reference levels of the radial stresses are assumed to be their values for the case of a rectilinear 
hollow formed between two sheets of an inextensible film of the thickness h under internal pressure 
p above ambient: 

_2pd 



( 14.35 ) 



On = 



Tch 



Stress Gi(0) can be calculated using ( 14.14 ) by the formula: 



( 14.36 ) 



a,(0) = 



liuzh 



14.7. Results of calculation 

Calculations were made by a program developed in the MathCAD-2(X)0 environment. The 
calculation used the following source data: 

d=lOcm,h = 0.01 cm, E = 36000 kg^cm^ v = 0.39, p =1 kg^cm^ 

Values of all dimensionless parameters listed in Section 14.6 were obtained for curvatures of the 
hoUow described by parameter r| in the range -10 < r| <10. 

11 problems were solved for the following values of parameter r|: 10"^ 10"° ^; 10"^ ^; 10"^"*; 10"^ ^; 1; 
10«-^; 10^^ 10'^ 10^^ 10. 

The results are listed in the table below. 

Table 3. Dimensionless parameters calculated using source data 

d=lOcm,h = 0.01 cm, E = 36000 kgf^cm^ v = 0.39, p =1 kgfi^cm^. 





e, 


02 


63 


sina 


Co 




C&=^& 




0.100 


1.068 


1.008 


0.992 


0.996 


0.100 


1.109 


0.925 


1.105 
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0.158 


1.098 


1.007 


0.982 


0.991 


0.159 


1.170 


0.889 


1.159 


0.251 


1.139 


1.005 


0.969 


0.980 


0.252 


1.269 


0.845 


1.244 


0.398 


1.192 


1.000 


0.952 


0.961 


0.397 


1.431 


0.797 


1.375 


0.631 


1.255 


0.994 


0.931 


0.930 


0.621 


1.693 


0.750 


1.575 


1.000 


1.321 


0.984 


0.908 


0.890 


0.967 


2.117 


0.708 


1.883 


1.585 


1.384 


0.973 


0.886 


0.845 


1.496 


2.795 


0.675 


2.361 


2.512 


1.437 


0.962 


0.867 


0.803 


2.302 


3.866 


0.652 


3.106 


3.981 


1.479 


0.953 


0.853 


0.773 


3.512 


5.536 


0.636 


4.279 


6.310 


1.510 


0.947 


0.843 


0.757 


5.298 


8.104 


0.625 


6.133 


10.000 


1.532 


0.943 


0.837 


0.756 


7.850 


11.998 


0.619 


9.073 



Shapes of the hollow's meridian for minimum, medium and maximum values of the curvature 
parameter r| are shown on Fig. 14.4, Fig. 14.5 and Fig. 14.6. The diagrams outhne the fold 
formation zone by a blue dotted line. As one can see from the diagrams, almost the whole shell 
except for a narrow region near the central glued area is covered by radial folds and is free from 
circular stresses. 
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Fig. 14.4. Shape of the meridian of the hollow at r| = 0.1 
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Fig. 14.5. Shape of the meridian of the hollow at T| = 1 
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Fig. 14.7. Dimensionless coefficients of shape 0i, 82 and 83 




Relative curvature 



Fig. 14.8. Sine of the peeling angle a of the seam of the central glued area 
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Fig. 14.9. Coefficients of concentration of: 
1 - radial stress in the hollow region near the boundary of the glued area; 
2 - linear peeling stress at the seam of the central glued area; 
3 - radial stress in interior points of the central glued area 
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Relative curvature 



Fig. 14.10. Coefficient of decrease of the radial stress and linear peeling stress 
in the glued area along the outer contour 

15. BIPEMODICAL MODELS OF AIRBAGS WITH ROUND GLUED AREAS 

This chapter features the results of the numerical computations of the stress and strain of endlessly 
stretching biperiodical models of airbags with round glued areas. 

The allotment of the linear peeling stress along the glued contour is investigated. 
Two types of models are investigated. In the first type model the centers of the round glued areas are 
located in the nodal points of a triangular grid, and in the second type model they are located in the 
nodal points of a square grid. 

The series of problems are resolved in which the glued area radius relations are varied to the grid 
spacing within the range 0.1, ..,0.3. 

The problem turned out to be unexpectedly difficult. There were observed, besides strong 
nonlinearity, the characteristics of the differential equations system degeneracy in the initial stage of 
the equations system deformation and degeneracy in the folds-forming stage (due to the degeneracy 
of the material itself). It was discovered at the last moment that the film deformable surface self- 
intersection is observed in the obtained solutions. This circumstance put the problem to the category 
of contact ones. 

As the numerical model is extremely unstable, there were taken some measures to stabilize it. The 
equibiaxial preliminary expansion of the airbag field was used as one of the stabilizing factors, 
which was gradually diminished in the process of the problem solution. As the other stabilizing 
factor an additional shell of special material with vanishingly small modulus of elasticity and very 
high coefficient of the thermal strain was used. 

15.1. The main calculated prerequisites 

The finite-element method, realized in multiphysical program system ANSYS, version 5.5.1. was 
used for computations. As the basic finite element the triangular membrane element SHELL41 is 
used which takes into account big displacements and deformations correctly and takes only tension 
stress. While working under conditions of negative deformations the element loses its response force 
in one or two directions imitating folds forming. At the same time the element becomes inwardly 
variable. 
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For this inward variability prevention, v^hich leads to the instability of the calculating process a 
workaround is used. It means that simultaneously with the basic shell of the finite elements the 
second shell of material with special characteristics is adjusted. This material has vanishingly small 
modulus of elasticity and very high coefficient of the thermal strain. The negative temperature is 
assigned for the elements of the second shell with the expectation of the thermal strain to be equal - 
1.0. An unstrained member of this material under the temperature influence shrinks to the zero size. 
At the same time low modulus of elasticity lets this element to be deformed almost without strain. 
Thus, the conglomerate of two elements - membrane and special - never undergoes compression 
stress; it is not intemally changeable and practically doesn't differ on its characteristics from the 
element SHELL41, which imitates folds forming. 

In order to stabilize the system on the initial stage of deformation when it is nearly flat as well as in 
the process of its skeleton reconstruction some efforts are made to the system, which causes 
proportional biaxial pretension. On the last stage of problem solution these efforts gradually reduce 
and in an ideal case should be reduced to zero. 

However, the moment comes when the shell heavily bulges and some parts of it exceed the plane of 
synmietry. Such deformation conflicts with the physical meaning as it corresponds to the strained 
surface with self-intersections. It results in the necessity to supplement the problem with the 
restricted inequations or to work in the rigid contact surface to the system, which coincides with the 
plane of symmetry. Such problem complication heavily increases (till hours) its solution time and 
brings additional elements of instability to the process of solution. 

In order to get the final problem solution we had to tum to rather rarefied finite elements system. 

15.2. Model with glued areas in nodes of a triangular grid 

The model being investigated results from two infinitely extending film sheets with die centres in 
nodes of a triangular grid being glued in round areas. The scheme of the work part for this model is 
shown on Fig. 15.1. 

Grid spacing a, round glued area radius r. Into the hollow between sheets gas is pumped under the 
pressure p. 




Fig. 15.1. The scheme of the work part for airbag with biperiodical system of round glued areas in 

the nodes of the triangular grid. 

Calculations were made with the following source data: 
• Film sheet sickness /i= 0.01 cm; 
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• Modulus of elasticity 36000 kg^cm^; 

• Poisson ratio v=0.39; 

• Grid spacing a= 10 cm; 

• Gas pressure above ambient p = 1 atm. 
Round glued area radius r was varied. 

The problems are solved for radii r: 1.0; 1.5; 2.0; 2.5; 3.0. 

Concidering symmetry characteristics present in this problem, it is enough to make a finite-element 
method solution only for unit cell, which has right triangle form shown on Fig. 15.1 and more 
detailed on Fig. 15.2. 




2 



Fig. 15.2. Unit cell scheme . 
1 - immovable planes of symmetry, 2 - movable plane of symmetry. 

On each side of this triangle mirror synunetry conditions are observed. Concidering missing fixation 
at infinity, one of the lines (or rather planes) of synmietry must be movable. 



Dummy finite elements of vanishingly small rigidity, which serve for gas pressure load transfer to 
the movable plane of symmetry, are located on the same plane. 

Unit cell splitting to finite elements is shown on Fig. 15.3. 




Fig. 15.3. Unit cell splitting to finite elements. Sizes are in millimeters. 

Airbags calculations with round glued areas of different radii are made with this splitting. Unit cells 
strained schemes are given on Fig. 15.5 - Fig. 15.9. Folds, connecting nearest round glued areas, are 
formed on the surface of airbag at small glued area radii . It can be seen on big finite elements 
contractions. 




Fig. 15.4. Typical sizes of airbag deformation. 

As the gas is pumped, the airbag inflates and its linear dimensions contract. 

Let the airbag thickness be indicated as d, and the size, which characterizes linear 
dimensions contraction as m. Geometrical meaning of these sizes and their connection with unit cell 
movement is evident from Fig. 15.4. 

Airbag thickness ratio dia and its contraction ula depending on glued area radius 

(15.1) p = -. 

a 

are listed in the Table 4. 

Table 4. Aiibag thickness ratio and contraction for the triangular grid case 





p = 0.10 


p = 0.15 


p = 0.20 


p = 0.25 


p = 0.30 


Thickness ratio dIa 


0.4323 


0.3635 


0.2947 


0.2285 


0.1682 


Contraction ula 


0.1026 


0.0798 


0.0572 


0.0372 


0.0209 
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Fig. 15.7. Unit cells strained scheme at glued area radius 2 cm. 
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Fig. 15.8. Unit cells strained scheme at glued area radius 2.5 cm. 




Fig. 15.9. Unit cells strained scheme at glued area radius 3 cm. 
Table 5. Linear peeling stress R45 for the triangular grid case in kgfi^cm. 





p = 0.10 


p = 0.15 


p = 0.20 


p = 0.25 


p = 0.30 


a= 0 


10.082 


7.557 


4.846 


3.362 


2.628 


a= 5 


2.576 


1.782 


1.308 


1.581 


1.374 


a= 10 


10.436 


5.386 


4.439 


3.124 


2.557 


a= 15 


12.747 


8.443 


5.586 


4.132 


3.179 


a = 20 


14.397 


9.442 


6.724 


4.879 


3.553 


a = 25 


15.086 


10.070 


7.200 


5.181 


3.738 


a = 30 


11.792 


8.508 


6.509 


4.927 


3.113 



Linear peeling stress is tabulated depending on the angle a (see Fig. 15.2) and relative glued area 
radius p. 

The same data are depicted graphically on Fig. 15.10. 
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Fig. 15.10. Linear peeling stress R^o for the triangular grid case. 

Linear peeling stress concentration coefficients k, which had been obtained by numerical method, 
are given in Table 6 . As a normalizing divisor an averaged hnear peeling stress is taken, which can 
be calculated by the following formula: 

(15.2) ^03„, W3a-tan30°_^^.Y 
2iir[ 2 ) 



Table 6.Linear peeling stress concentxation coefficients ^30 
for the triangular grid case . 





p = 0.10 


p = 0.15 


p = 0.20 


p = 0.25 


p = 0.30 


a= 0 


0.759 


0.896 


0.822 


0.789 


0.849 


a= 5 


0.194 


0.211 


0.222 


0.371 


0.444 


a=10 


0.786 


0.638 


0.753 


0.733 


0.826 


a= 15 


0.96 


1.001 


0.948 


0.969 


1.027 


a = 20 


1.084 


1.119 


1.141 


1.144 


1.148 


a = 25 


1.136 


1.193 


1.222 


1.215 


1.208 


a = 30 


0.888 


1.008 


1.105 


1.156 


1.006 



In order to get a real value of the Unear peeling stress R^o, concentration coefficient from Table 6 
must be multiplied by the normalizing divisor calculated by formula ( 15.2 ): 

( 15.3 ) R^Q = k^oR^Q. 



15.3. Model with glued areas in nodes of a square grid 

Calculations for the model with round glued areas in nodes of a square grid are made similarly to the 
previous case. 

The scheme of the work part for this model is shown on Fig. 15.1 1. 
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Fig. 15.1 1. The scheme of the work part for airbag with biperiodical system of round glued areas in 

the nodes of a square grid. 

As in the previous case, grid spacing a, round glued area radius r. Into the hollow between sheets 
gas is pumped under the pressure p. 

Calculations were made with the following source data: 

• Film sheet sickness h- 0.01 cm; 

• Modulus of elasticity E= 36000 kg^cm^; 

• Poisson ratio v=039; 

• Grid spacing a= 10 cm; 

• Gas pressure above ambient p- \ atm. 
Round glued area radius r was varied. 

The problems are solved for radii r: 1.0; 1.5; 2.0; 2.5; 3.0 cm. 
The form and sizes of the unit cell are shown on Fig. 15.12. 
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Fig. 15.12. Unit cell scheme . 
1 - immovable planes of synmietry, 2 - movable plane of symmetry. 

Unit cell splitting to finite elements is accepted the same for all problems. 
It is shown on Fig. 15.13. 




Fig. 15.13. Unit cell splitting to finite elements. Sizes are in millimeters. 
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Airbags calculations with round glued areas of different radii are made with this splitting. Unit cells 
strained schemes are given on Fig. 15.14 - Fig. 15.18. Deep folds, connecting nearest round glued 
areas, are formed on the surface of airbag at small glued area radii . It can be seen on the pictures. 

. At big glued area radii the surface of the airbag remains practically smooth. 




Fig. 15.14. Unit cells strained scheme at glued area radius 1 cm. 




Fig. 15. 15. Unit cells strained scheme at glued area radius 1.5 cm. 
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Fig. 15.16. Unit cells strained scheme at glued area radius 2 cm. 




Fig. 15.17. 

Unit cells strained scheme at glued area radius 2.5 cm. 
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Fig. 15.18. Unit cells strained scheme at glued area radius 3 cm 
As the gas is pumped, the airbag inflates and its linear dimensions contract. 

These deformations are evaluated quantitatively by the relative airbag thickness dia and its 
contraction w/a, which is uniform in all directions (see Fig. 15.4). 



Table 5. Airbag thickness ratio and contraction 





p = 0.10 


p = 0.15 


p = 0.20 


p = 0.25 


p = 0.30 


Thickness ratio 


0.5566 


0.5027 


0.4253 


0.3486 


0.2750 


Contraction 


0.1178 


0.1054 


0.0798 


0,0565 


0.0360 



Table 7. Linear peeling stress /?45 for the square grid case in kgf/cm. 





p = 0.10 


p = 0.15 


p = 0.20 


p = 0.25 


p = 0.30 


a= 0 


12.381 


8.410 


5.617 


3.860 


2.688 


a= 5 


-0.360 


-0.992 


-0.153 


0.828 


1.180 


a= 10 


4.274 


2.764 


2.235 


2.130 


1.544 


a= 15 


6.493 


3.581 


2.822 


2.411 


1.751 


a = 20 


9.888 


5.835 


4.185 


2.885 


2.331 


a = 25 


13.438 


9.619 


6.602 


4.682 


3.703 


a = 30 


17.429 


12.283 


8.826 


6.485 


4.785 


a = 35 


18.942 


14.015 


9.871 


7.125 


5.343 


a = 40 


18.481 


13.574 


10.135 


7.541 


5.650 


a = 45 


16.591 


12.277 


9.244 


6.917 


5.511 



The same data are shown graphically on Fig. 15.19. 
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alpha, degrees 

r=1.0 

r=1.5 

r=2.0 

r=2.5 

r=3.0 



Fig. 15.19, Linear peeling stress R45 for the square grid case . 

Linear peeling stress concentration coefficients fcjs, » which had been obtained by numerical 
method are given in Table 8 . 

As a normalizing divisor an averaged linear peeling stress is taken which can be calculated by the 
following formula: 

(15.4) R%,^-P-(a^-nr^) 



Table 8. Linear peeling stress concentration coefficients fcts 
for the square grid case . 





p = 0.10 


p = 0.15 


p = 0.20 


p = 0.25 


p = 0.30 


a= 0 


0.803 


0.853 


0.807 


0.755 


0.706 


a= 5 


-0.023 


-0.101 


-0.022 


0.162 


0.310 


a= 10 


0.277 


0.280 


0.321 


0.416 


0.406 


a= 15 


0.421 


0.363 


0.406 


0.471 


0.460 


a = 20 


0.641 


0.592 


0.601 


0.564 


0.613 


a = 25 


0.872 


0.976 


0.949 


0.915 


0.973 


a = 30 


1.131 


1.246 


1.269 


1.268 


1.257 


a = 35 


1.229 


1.421 


1.419 


1.393 


1.404 


a = 40 


1.199 


1.377 


1.457 


1.474 


1.485 


a = 45 


1.076 


1.245 


1.329 


1.352 


1.448 



Page 161 



In order to get a real value of the linear peeling stress R45, concentration coefficient from Table 8 
must be multiplied by the normalizing divisor calculated by formula: ( 15.4 ): 

( 15.5 ) /?45 = /l45/?^45. 



It is necessary to note that the problem to be solved is ill-conditioned regarding its physical nature. 
Obviously, full system energy has many local minimums and iterative process finds one of them 
occasionally. That explains the evident scatter of readings, given in tables and diagrams. 

16. AIRBAG WITH NARROW CHANNELS BETWEEN HOLLOWS 

V- 

This chapter discusses results of a numerical analysis of the stress and strain distribution in 
periodical models of airbags having narrow channels between their adjacent hollows. 

Principal properties of the shape of an inflated airbag and the distribution of the linear peeling stress 
along the glued contour have been studied. 

16.L Key assumptions of the analysis 

The analysis used a finite element method implemented in a multi-physical computational software 
ANSYS, version 5.5.1. The main finite element used was a triangular membrane element SHELL41 
that allows for large displacements and strains correctly and resists to tension only. 

16.2. A model of an airbag having channels between its adjacent hollows 

A plane schematic of a few cells of the model is shown in Fig. 16.1. As can be seen from the figure, 
the model in question is made of two infinitely extending strips of film 2(a-^b) wide glued in the 
shaded areas. The radius of the glued area is r, the width of the channel is 2z. The glued area and the 
whole system possess a translation symmetry with spacing 2c. 




Fig. 16.1. A schematic of cells of the periodic model of an airbag with air channels 

between adjacent hollows. 

A gas under pressure p is pumped to the hollow between the two film sheets. 

The following source data were used for the analysis: 

• Fihn sheet thickness h= 0.01 cm; 

• Elasticity modulus E= 36000 kg/cm^; 

• Poisson ratio v=0.39; 

• Gauge pressure of the gas p = 1 atmosphere. 
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• Fixed geometric sizes: a= 5 cm; 6=0 cm; r=l cm. 
Sizes c and z were subject to variation. 

The problems were solved with c equal to 3; 5; 7; 9 cm at z = 1 cm and with z equal to 0.5; 1.0; 1.5; 
2.0 cm at c = 5 cm. 

The value of b was assumed to be zero as one that had little effect on the result of the analysis. 

Due to the properties of synmietry possessed by the shape, it suffices to perform the finite element 
analysis only for an elementary cell of the rectangular shape shown in Fig. 16. land in Fig. 16.2 
more detail. 
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Fig. 16.2. A schematic of an elementary cell 
(1) fixed planes of synunetry, (2) moving plane of symmetry 



Conditions of bilateral symmetry hold at three sides of this rectangle. As there is no fixation in the 
infinity, one of the symmetry planes must be capable of moving. The arrangement of the synmietry 
planes is shown in Fig. 16.2. 

The division of the elementary cell into finite elements is shown in Fig. 16.3. 
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Fig. 16.3. A division of an elementary cell into finite elements, 
a = 5 cm; c = 9 cm; r = 1 cm; z = 1 cm. 
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16.3. Dependence of sizes of the airbag and the peeling stress at the boundary of the 
glued area on the length of the hollow 

Given below are results of solution of a series of problems where the length of the hollow was a 
variable value. Fig. 16.4 to Fig. 16.7 show deformed shapes of the elementary cell and give principal 
source data used in the analysis. 




Fig. 16.4. A deformed shape of the elementary cell 
at a = 5 cm; c = 3 cm; r = 1 cm; z = 1 cm. 




Fig. 16.5. A deformed shape of the elementary cell 
at a = 5 cm; c = 5 cm; r = 1 cm; z = 1 cm. 
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Fig. 16.6. A deformed shape of the elementary cell 
at a = 5 cm; c = 7 cm; r = 1 cm; z = 1 cm. 



fl2 




Fig. 16.7. A deformed shape of the elementary cell 
at a = 5 cm; c = 9 cm; r = 1 cm; z = 1 cm. 

These schematics show that folds are formed in the areas of the channels. The folds appear also in 
the glued region, which results in this region narrowing in the direction of the hollow's comer. 

The geometric dimensions of the protective airbag in its inflated state have been studied: thickness 
d, the channel's height/, contractions Ux and Uy, These are shown in Fig, 16.7. 

Values of these dimensions are given in Table 4. 



Table 9. Geometric dimensions of the protective airbag in its inflated state 





c = 3 cm 


c = 5 cm 


c = 7 cm 


c = 9 cm 


Airbag thickness d 


2.9282 


4.6068 


5.4414 


5.9452 


Channel height / 


1.0790 


1.1826 


1.2652 


1.2082 
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Contraction in x direction 


0.6579 


1.0583 


1.1843 


1.1921 


Contiraction in y direction 


0.6530 


1.1295 


1.4567 


1.6153 



Fig. 16.8 shows dependencies of the airbag thickness d and the channel height /on the c dimension. 




0.2 



^3456789 ^3456789 
c, cm c, cm 

Fig. 16.8. Dependence of the airbag thickness d and the channel height/on the c size 

The thickness of the airbag grows as the hollow length increases, tending to the value 20/n = 6.366. 
The channel height tends to the value 4/k = 1.273. Some instability in tfie behavior of this quantity 
can be noticed. The explanation for this is an instability of the fold formation. 

Fig. 16.9 shows diagrams of relative contractions of the airbag in the directions x and y. 




0.05 

^3456789 
c 

X direction 

y direction 

Fig. 16.9. Relative contractions of the airbag 
The maximum peeling stress can be detected in the semi-circular glued area. This is tabulated vs. 
angle a (see Fig. 16.2) and size c in Table 10. 

Table 10. Linear peeling stress in the semi-circular glued area, kgf^cm. 





c = 3 cm 


c = 5 cm 


c = 7 cm 


c = 9cm 


a= 0 


0.4585 


0.4481 


0.5788 


0.6602 


a= 10 


0.6166 


0.6642 


0.5736 


0.6447 


a = 20 


0.5496 


0.5633 


0.6269 


0.4201 


a = 30 


0.5026 


0.6006 


0.3301 


0.6149 
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a = 40 


0.4797 


0.8533 


0.2201 


0.2688 


a = 50 


2.3324 


2.9152 


2.9536 


3.6338 


a=60 


6.9014 


7.2006 


10.1490 


8.9788 


a = 70 


6.7891 


11.7564 


9.5840 


8.1324 


a=80 


8.0871 


8.7404 


10.6711 


14.9662 


a = 90 


1.2990 


2.0780 


2.0958 


1.6554 



The same data are shown graphically in Fig. 16.10. 




alpha, degrees 

c = 3 

c = 5 

c = 7 

c=9 

Fig. 16.10. Linear peeling stress in the semi-circular glued area 

As can be seen on the diagrams, at angles 60® < a < 80° a noticeable stress concentration appears. 
The spread of the values can be explained by the contiguity of the folded region. 

16.4. Dependence of the airbag dimensions and the peeling stress at the boundary of 
the glued area on the channel width 

Given below are results of solution of a series of problems where the variable quantity is the channel 
width z. Fig. 16.5 to Fig. 16.14 show the deformed shapes of the elementary cell and give principal 
source data used in the analysis. 




Fig. 16.1 1 . A deformed shape of the elementary cell 
at a = 5 cm; c = 5 cm; r = 1 cm; z = 0.5 cm. 
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Fig. 16.12. A deformed shape of the elementary cell 
at a = 5 cm; c = 5 cm; r = 1 cm; z=l cm. 




Fig. 16.13. A deformed shape of the elementary cell 
at a = 5 cm; c = 5 cm; r = 1 cm; z = 1.5 cm. 




Fig. 16.14. A deformed shape of the elementary cell 
at a = 5 cm; c = 5 cm; r = 1 cm; z = 2 cm. 

The geometric dimensions of the protective airbag in its inflated state have been studied: thickness 
dj the channel's height /, contractions Ux and Uy. These are shown in Fig. 16.7. 
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Values of these geometric dimensions are given in Table 1 1 . 



Table 1 1 . Geometric dimensions of the protective airbag in its inflated state 





7 = 0 S cm 


7=10 cm 


7 = 1 S f*m 

^ M. Will' 


7 = 20 cm 

^~ ^*\J will 


Airbag thickness d 


4.5284 


4.6068 


4.7850 


4.9924 


Channel height / 


0.5882 


1.1826 


1.3356 


2.3680 


Contraction in x direction 


1.1514 


1.0583 


1.0705 


1.0822 


Contraction in y direction 


1.1232 


1.1295 


1.1711 


1.1988 



Fig. 16.15 shows the dependence of the airbag thickness d and the chaimel height /on size z. 




Fig. 16.15. Dependence of the airbag thickness d and the height of the channel/on size z 

The thickness of the aiibag increases a bit as the channel becomes wider. The height of the channel 
grows as its width increases. In the diagram of Fig. 16.15 one can see a non-smooth variation of the 
value of/. The explanation for this is that it is the channel's vicinity where the folds arise, and the 
central point of the chaimel can be sometimes at the top of a fold and sometimes at the bottom of a 
concave. 

Fig. 16.16 shows diagrams of the relative contraction of the aiibag in the direction of x and y. 

0.35 I 

0.3 

0.25 ' ,,J.. 

0.2 ' ' 

0.15 

0.1 

0.05 1 

%,5 1 1.5 2 

z, cm 

X direction 

y direction 

Fig. 16.16. Relative contraction of the airbag 
These contraction values depend litde on the width of the channel. 

The maximum Unear peeling stress can be observed in the semi-circular glued area. This is tabulated 
vs. angle a (see Fig. 16.2) and size z in Table 12. 

Table 12. Linear peeling stress in the semi-circular glued area, in kgfi^cm. 
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z = 0.5 cm 


z = 1 cm 


z = 1.5 cm 


z = 2 cm 


a= 0 


0.3163 


0.4481 


0.2327 


1.4292 


a= 10 


0.3840 


0.6642 


0.4252 


2.1043 


a = 20 


0.3100 


0.5635 


0.6430 


2.9817 


a = 30 


0.1501 


0.6006 


0.9163 


3.1759 


a = 40 


0.0791 


0.8533 


1.1897 


5.1840 


a =50 


1.0590 


2.9146 


8.1668 


8.3181 


a = 60 


4.0550 


7.2006 


7.6596 


7.7622 


a=70 


6.9547 


11.7564 


8.6510 


8.1868 


a=80 


12.1089 


8.7404 


13.5679 


6.6607 


a = 90 


1.7862 


2.0780 


0.3443 


1.0787 



The same data are presented graphically in Fig. 16.17. 




Fig. 16.17. Linear peeling stress in the semi-circular glued area 

As one can see on the diagrams, at angles 50° < a < 80° a substantial stress concentration appears. 
The spread of the values is even greater than that in Fig. 16.17. It can be explained by the contiguity 

of the folded area. 

The analysis enables us to have only a qualitative idea of the magnitude of the peeling stress. This 
stress reaches the values of 8 to 15 kgfi^cm^. 

17. AIRBAG HAVING A LONG GAS-FEEDING HOLLOW 

This chapter presents results of a numerical analysis of the stress/strain distribution in periodical 
models of airbags having long hollows for the purpose of gas feeding. 

The objects of investigation include qualitative properties of the shape of an inflated airbag and a 
distribution of the linear peeling stress in the seam of the elliptic adhesive joint. The effect of stress 
concentration in an area of a maximum curvature of the ellipse has been studied in detail. 
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17.1. Analytical tools employed 

The analysis was performed using the finite element method as implemented in the multi-physical 

software system ANSYS, version 5.5.1. The principal finite element used was a triangular 
membrane element SHELL41 that allowed for large displacements/strains and resisted to tension 
only. 

17.2. Models of adjoining and intersecting hollows 

The airbags implement a popular design solution in which a long hollow commonly used to feed the 
gas is adjoined by a few shorter hollows perpendicularly or at another angle. There are areas of 
stress concentration in places where the hollows join one another. 




_d 
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A 



f 



d 



Fig. 17.1. An airbag model with a long gas-feeding hollow 

A simpUfied model of this design and its elementary cell are shown in Fig. 17.1. 
Same investigations can be performed with a bi-periodical model shown in Fig. 17.2. In this case the 
elementary cell is smaller in size and can be covered by a denser mesh of finite elements. This 
model possesses an advantageous difference fi'om the previous one: in its deformed shape it 
maintains a translation symmetry [ 19 ] while that model loses this kind of symmetry when 
deformed. The quantities in question are affected by the model type only a littie bit. 

The model shown in Fig. 17.2 was chosen for the serial analysis. 
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Fig. 17.2, An airbag model with a long gas-feeding hollow and a smaller elementary cell 
173. Description of the series of problems solved 

The design model includes an eUiptic adhesive joint. The ratio of the ellipse's semi-axes influences 
the distribution of the peeUng stress at the boundary of the adhesive joint. The more extended the 
ellipse is, the greater is the stress concentration in the place of a maximum curvature of its boundary. 
It is exactly this phenomenon that this chapter deals with. A series of problems was solved at a 
constant size of the semi-major axis (3 cm) and a variable size of its semi-minor axis e. The analysis 
was performed at ^ = 1.00 cm; 1.25 cm; 1.50 cm; 1.75 cm; 2.00 cm. The other parameters are the 
same in all problems: 

the film sheet thickness h- 0.01 cm; 
elasticity modulus E- 3600 kg/cm^; 
Poisson ratio v=0.39; 
half-width of the main hollow a= 5 cm; 
half-width of the distributive hollow c= 3 cm; 
half-width of the channel z- 3 cm; 
semi-major axis of the eUipse 3 cm; 
length of the elementary cell 15 cm; 
gas pressure above ambient = 1 atmosphere. 
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All these quantities are shown in a detailed pattern drawing of the elementary cell presented in Fig. 
17.3. 

The choice of a lower elasticity modulus E- 3600 kg/cm^ should be emphasized. Numerous 
attempts undertaken by the author to obtain a stable solution for a stiff material with its elasticity 
modulus E = 36000 kg/cm^ had no success. In all cases folds were formed causing a saw-tooth 
distribution of the peeling stress along the contour of the adhesive joint. The shape of the folds was 
determined apparentiy by the method of meshing into finite elements rather than the physical sense 
of the problem. Having understood this, the author lowered ten times the stiffness of the material 
thus obtaining smootii shapes of the bag in its inflated state and smootii curves of the peeling stress 
distribution along the adhesive joint. It made possible to analyze the results of the calculation. There 
should be no big error caused by the lower stiffness of the material because in both cases the peeling 
stress is integrally balanced by nearly the same gas pressure. 




/=15 



/=20 



E 

o 
in 



15 cm 



Fig. 17.3. A schematic of the elementary cell. The dashed line shows the boundary of the adhesive 
joint that resists to the peeling stress. Also, the principle of the point numbering at the elliptic joint 

boundary is shown 

17.4. Results of the solution 

All five problems have been solved successfully. 

Fig. 17.4 to Fig. 17.8 show deformed shapes of the elementary cell at different values of e. 
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Fig. 17.4. A deformed shape of the elementary cell at e = 1 .00 cm 





Fig. 17.7. A deformed shape of the elementary cell ate= 1.75 cm 
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Fig. 17.8. A defonned shape of the elementary cell at ^ = 2.00 cm 




Fig. 17.9. A distribution of the peehng stress along the boundary of the elliptic adhesive joint 

in more extended ellipses 
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Fig. 17.10. A distribution of the peeling stress along the boundary of the elliptic adhesive joint 

in less extended ellipses 

Also, the peeling stress at the boundary of the elliptic adhesive joint has been determined by this 
analysis. The boundary was divided into 20 segments of an equal length. End points of the segments 
are numbered from 0 to 20 as shown in Fig. 17.3. 

The peeling stress is calculated in these points and represented as plots in Fig.9 and Fig. 10. 

The same diagrams show the level of the peeling stress for two idealized cases as horizontal lines. 
The first one of those (h = 10 cm) conforms to the case of an infinite cylindrical hollow glued out of 
two strips 10 cm wide (the main hollow), the second case (h = 6 cm) corresponds to the case of an 
infinite cylindrical hollow glued out of two strips 6 cm wide (the gas distributing hollow). 
As the plots demonstrate, the level of the stress in the intersecting hollows is always higher than that 
in the idealized case of non-intersecting hollows of the same size. The more extended the elliptic 
joint, the higher the peeHng stress at its boundary. As the length of the semi-minor axis of the ellipse 
increases, the peeling stress at the vertex of the ellipse falls. This dependence is shown as a plot in 
Fig. 17.11. 




1 1.25 1.5 1.75 2 

small bal£axe of ellipse, cm 



Fig. 17.1 1. Peeling stress at the ellipse's vertex (/ = 10) 
V. its semi-minor axis 
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At the beginning {e < 1.25 cm) the peeling stress at the vertex of the ellipse exceeds that in nearby 
points. Then (e > 1.50 cm), the peeling stress at the vertex becomes less than that in the nearby 
points. Apparently, the rational size of the semi-minor axis of the eUipse in the presented series of 
problems should be 1.40 cm or so, when the stress at the ellipse's vertex becomes level with that in 
points located near the vertex. 
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